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ABSTRACT 


This thesis investigates the periodic autocorrelation function (PACT) and periodic 
ambiguity function (PAF) for orthogonal continuous waveform (CW) modulations used 
in netted low probability of intercept (LPI) radar. Three orthogonal polyphase sequences 
and one frequency coding sequence are examined and their PACF and PAF 
characteristics are quantified. The Wigner-Ville distribution (WVD) and quadrature 
mirror filter bank (QMFB) time-frequency signal processing techniques and the 
cyclostationary bi-frequency technique (often used in non-cooperative intercept 
receivers) are used to detect the orthogonal CW signals and extract their parameters. The 
results shows that a combination of the techniques used were able to extract the basic 
signal parameters of bandwidth and code period from the polyphase waveforms and also 
the frequency hop slots and code length from the frequency coding sequence. The 
concept of using a swarm of unmanned aerial vehicles (UAV) is examined from the 
viewpoint of a coordinated group of netted intercept receivers in search of an LPI radar 
network. 
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I. INTRODUCTION 


A. ORTHOGONAL NETTED LPI RADAR SYSTEMS 

The term LPI is that property of radar theory that, because of its lower transmitted 
power, wide bandwidth, frequency variability, or other design characteristics, make it 
difficult to be detected by a non-cooperative passive intercept receiver. The principle of 
target detection by the radar prior to the target being able to detect the electromagnetic 
energy from the radar is the key to LPI. Radars that possess these qualities reduce their 
vulnerability to detection, location, and most importantly destruction. LPI radars are able 
to provide the operator with warning of incoming targets while remaining relatively 
invisible to the target. This distinct advantage is the impetus for the great interest in the 
ability to produce viable LPI radars and systems that can detect an adversarial LPI radar 
system. LPI radar is typically identified by one or more of the following characteristics 
[ 1 ]: 

• Low average power continuous wave (CW) radiation with a large time 
bandwidth product. 

• Radiated energy that is spread over a wide angular region, wide frequency 
band, and long time interval. 

• Ultra-low side lobe transmit antenna 

• Power management via radiation only as necessary to obtain target 
information, reduced levels to maintain situational awareness of target 
movement, and a reduction in power as a target gets closer to the radar. 

• Utilization of the properties of the atmosphere to attenuate a radar signal 
beyond a usable range. 

The networking of LPI radar systems will further reduce their vulnerability 
probability of detection by non-cooperative electronic support receivers (ESR). A 
networked LPI system can share target information in one of two ways. The first method 


1 



would simply share a common operating picture, passing target data between LPI radars 
to fuse the data and refine the target. The second method would rely on each radar 
system transmitting its own signal, but process echoes from the target from its own signal 
and also process echoes from the other LPI transmitted waveforms. This system could 
reduce its total transmitted power requirement as it would integrate the returns from the 
other LPI radars to maintain track on the target. Both cases present a difficult target for 
anti-radiation missiles. 

The study of LPI radars requires two perspectives to be effective. This concept is 
to ensure that the LPI radar provides an operator with target information while remaining 
below the threshold of the adversarial ESR ability to identify the LPI radar. The other 
perspective is that of the operator of the ESR, to identify and locate an adversary’s LPI 
radar. This thesis looks through both “lenses” first by incorporating advances in signal 
design to a networked LPI system and second by utilizing a networked swarm of UAV 
equipped with ESR capabilities that incorporate advanced signal processing to identify 
and locate the adversary’s LPI radar. 

B. PRINCIPAL CONTRIBUTIONS 

The first objective of this thesis is to examine several orthogonal LPI radar 
waveforms and their vulnerability to the detection by a non-cooperative ESR utilizing 
time frequency and bi-frequency detection methods. The second objective is to illustrate 
the concept of swarm theory as applied to a network of UAV equipped with the non- 
cooperative ESR to show the benefits of a stand-in sensor approach to LPI radar 
detection. 

The orthogonal LPI signals were developed in MATLAB [2] and based on code 
sequences [3-5] in a program titled “LPI Signal Toolbox” [6], which also includes other 
contributions by the authors of this thesis. The netted LPI radar signals generated include 
three orthogonal polyphase coded (OPS) CW waveforms and a discrete frequency coding 
(DEC) or frequency hopping waveform was simulated. 

Radar receiver performance was determined by examining the autocorrelation 

function (ACF), periodic autocorrelation function (PACF) and periodic ambiguity 
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function of the LPI waveforms. Non-cooperative receiver processing techniques 
including the Wigner-Ville Distribution, Quadrature Mirror Filter Bank, and 
Cyclostationary Spectral Analysis [1] were applied to study how the waveform 
parameters were revealed in the time-frequency and bi-frequency domains. The concept 
of swarm theory using a UAV as a stand-in sensor to locate and identify the netted LPI 
radar systems is studied. The simulation of the radar detecting UAV swarm in Java Code 
shows the effectiveness of networking a swarm of UAV to enhance detection capability 
of the netted LPI radar system. 

C. THESIS OUTLINE 

Chapter II provides a brief description of LPI receiver analysis techniques used 
and the mathematical principles that govern LPI radar performance. This description is 
necessary as these techniques quantify the performance of each signal type and allow for 
a direct comparison of their characteristics. 

Chapter III presents the orthogonal netted radar waveforms, a description of their 
characteristics and the MATLAB® [2] tool that was used to create the simulated 
waveform. Receiver characteristics of each waveform are described and compared. 

Chapter IV shows the results of the time-frequency and bi-frequency detection 
techniques used in the ESR. 

Chapter V introduces swarm theory and the advantages a UAV Swarm has in 
detection and location of networked LPI radars. 

Chapter VI will contain a brief discussion of the JAVA program code for 
simulating the UAV swarm, instructions for use of the program, and a battlefield 
application. 

Chapter VII will summarize this thesis and provide recommendations for future 
research. 
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II. NETTED RADAR RECEIVER ANALYSIS 


In this chapter, the signal-to-noise ratio (SNR), the periodie ambiguity funetion 
(PAF), and the periodic autocorrelation funetion (PACF) of a system of netted LPI radars 
is discussed. The cross correlation funetion (CCF) used to evaluate the orthogonality of 
the netted radar waveforms is also diseussed. 


A. SIGNAL-TO-NOISE RATIO 


The power reeeived by any single radar is: 

_ P,G,G^X^(7 


( 1 ) 


where, is the reeeived power, P, is the transmitted power, G, is the transmit antenna 


gain, is the receive antenna gain, A is the transmitted signal wavelength, cr is the 


radar eross seetion of the target, is the range from the radar transmit antenna to the 


target, R^ is the range from the target to the reeeive antenna (the two ranges are equal for 

a monostatic radar, but in netted radar theory, the transmitter and receiver are often not 
co-loeated). The received power also faetors into the SNR ealculation. The SNR is often 
used to determine the threshold at which the radar or netted radars will be able to detect 
an incoming target [9]. 


In most networked radar systems, the nodes in a network eonsist of transmitters 
and reeeivers that have been engineered to be very similar if not eompletely alike. Given 
the likelihood of eaeh node having very similar charaeteristies, the SNR equation is 
altered to isolate the immutable parameters from those that can be ehanged by [9]: 

1 


SNR = ^ Y 


( 2 ) 


where t is the target integration time in seconds, k is Boltzmann’s constant, is the 


system temperature, is the noise factor of the system, L is the term that ineorporates 
any system losses in the signal processing, m and n are summation indiees that aeeount 
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for all the transmit and receive antennas in the radar network. From this new form of the 
SNR equation, it is evident that range to the target greatly affects the received SNR. 
Most radar designs have the ability to adjust transmitted power as needed to acquire their 
targets at maximum detection range. However, the main premise of LPI radar is to 
maintain the transmit power at a very low level (on the order of watts or milliwatts) to 
avoid detection by an adversary intercept receiver. Hence, a target’s range is related 
indirectly to the receive SNR [18]. 


The ability of netted radars to collectively assess the position of the target due to 
the multiple received signals is the main advantage over single radars. Given N co¬ 
operative radars in a battle space, there are N possible radar returns that any individual 
radar in the network can process. Assuming an isotropic radiator as the target and the 
noise as the primary limiting factor in network performance, the netted SNR Equation (2) 
normalized to the case of a single radar ((^SNR)^ ) yields [18]: 


(SNR)ff 1 

(SNR), 


(3) 


This result takes into account the removal of all the common terms in the equation. 


Assuming thati?,, = Kj ^ Equation (3) is proportional to the square of the number of 


individual radars in the network, or [18]: 


(SNR), 


(4) 


B. THE AMBIGUITY FUNCTION 


The ambiguity function (AF) of all the radars in a network can be weighted in 
order to achieve a more accurate target position. For a simplistic view of this concept, it 
is assumed that: the noise levels at the receivers are all the same, the transmitted signals 
and matched fdters are all the same signifying a coherent system, the returns from the 
target are not aspect dependent and therefore the target can be treated as an isotropic 
radiator. With these assumptions, the weighting factors for all N radars in the network 
can be calculated using the received power equation and comparing that to the maximum 


6 



received power by any of the radars in the network. The weighting factor is determined 
by comparing the received power of one radar Equation (1) with the maximum received 
power of any radar in the network using 

w.=—^ - (5) 

Max{P^) 

with w. as the weighting factor for the ith radar, and Max[P^) is the maximum received 
power of any receiver in the network [9]. 


The AF is an integral part of the radar’s ability to process target detections using a 
low transmitted power level. The theory behind the AF is the use of a copy of the 
transmitted signal to identify anomalies (alterations in the transmitted signal) that 
correspond to potential target echoes. A matched fdter is used to process the incoming 
received signal. The matched fdter is directly correlated to the transmitted signal. When 
the received signal is processed through the matched fdter, the differences in the Doppler 
shift of the received signal correspond to potential target detections. A positive v 
implies a target is moving towards the radar receiver. A positive range offset t means 
the target is further from the radar than the reference value ofr = 0. The ambiguity 
function should be viewed as the envelope of the matched fdter output when the received 
signal is the input (even this signal is just the transmitted signal with any Doppler shifts 
that correspond to potential targets). The equation below represents the AF [8]: 




t-UJ 

I u{t)u {t + dt 


( 6 ) 


The combination of the weighting factor from Equation (5) and the individual AF results 
from the individual radars Equation (6) can be summed to provide an overall netted 
ambiguity for the target or targets of interest [9]: 


^netted 




i=l 


( 7 ) 


There are three properties of the AF that are the basis of the analyses discussed 
later in this paper. The first of these defines the origin as the location for the maximum 
of the AF. This is obtained by normalizing the received signal energy, which in turn 
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yields a value of unity at the origin| 2 '(r,v)| < |j(0,0)| = 1 [10]. The seeond property 

refers to the conservation of the energy of a signal. Specifically, the total volume under 
the square of the normalized ambiguity surface is equal to unity and is independent of the 

particular signal’s waveform as shown by [10]: 

00 00 

dTdv = \ ( 8 ) 

—00 —00 

The third property describes the symmetry of the AF with respect to the 
origin|j(-r,-v)| = |j(r,v)| [10]. This allows for the focus of analytical efforts on the 

two adjacent quadrants of the AF. The analyses and plots will only contain quadrants 1 
and 2, which will show the positive Doppler shift values of the received signal. 

C. THE PERIODIC AMBIGUITY EUNCTION 

The AF is defined for signals of a finite duration, so it can be surmised that the 
periodic ambiguity function (PAF) uses multiple periods of the AF and is therefore 
defined for infinite duration signals. The LPI signals discussed in this paper are CW 
waveforms, so the PAF will be used to analyze these signals of infinite duration in the 
same way the AF is used to analyze the pulsed waveforms of finite duration [1]. The 
PAF is shown below as the normalized response of the matched received that correlates 
an infinite length of received signal with a reference signal of finite duration. The 
equation for a single reference code PAF would be a slight modification of the AF 
Equation (6) [10]: 

Xr(T.v)\=^\u(t)u*(t+Ty‘^--'"dt (9) 

0 

where u{t^ is the received CW signal of infinite duration and is the received signal 
code period. The reference signal of A’ code periods have a PAF as shown below [14]: 

Z„ (r,v) = — f u(t)u*(t + T)e‘^’"dl (10) 

'■ N 1 


8 



This is N’ copies of the reference signal being integrated over the length of the received 
signal. Using an important relationship between Equation (9) and Equation (10) [14]: 


kivT(uv)| = |jr(r,v)| 


sin(^N'7rvT^) 

N'sm(^7rvT^) 


( 11 ) 


This suggests a simplified process to obtain the N’ period PAF from the single period 
PAF by multiplying the single period PAF by [14]: 


A'sin(;Tv7].) 


( 12 ) 


D. THE AUTOCORRELATION EUNCTION 


The autocorrelation function (ACF) is generally used for signals of finite energy. 
The incoming signal is detected by the receiver and sent to its matched filter. The 
received signal, which is the same as the transmitted signal except for any target 
interactions, is then shifted along the time axis by varying time shifts. During the shifting 
process of the received signal, the transmitted signal is correlated with the time shifted 
received signal to identify any anomalies present. These differences of the transmitted 
and received signals correlated with one another over the time interval defined by the 
matched filter can be attributed to targets detected by the radar. 


\R(t)\ = 


I u(t)u*(t + T)dt 


= kr(u0)| 


(13) 


As shown above, there is a unique relationship between the two-dimensional AF with the 
Doppler shift equal to zero and the ACF. A Doppler shift of zero equates to the 
equivalent range window of the matched filter in the receiver. This is only a single 
instance of the matched filter having the same range window as the ACF, but it illustrates 
the relationship between the two functions and their properties [10]. 


E. THE PERIODIC AUTOCORRELATION FUNCTION 


The PACF is related to the PAF in the same way the ACF is related to the AF. 
The CW waveforms that are typical of LPl signals are often modulated by a periodic 
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function such as a phase code sequence or a linear frequency ramp function. The 
equation for the PACT shows its reliance on subcode duration delay r, and the 

complex envelope m (t): 

c i=l 

The PACT is also inversely proportional to the number of subcodes present in the 

phase code sequence, which further annotates the distribution of the signal’s energy. The 
PACF is also the zero Doppler cut of the PAF as shown by the 

relationship (r,0) =i?(rt^). LPI waveforms can sometimes yield a perfect PACF 
where the function reaches unity at t = 0,T^,2T^,... and is zero otherwise [11]. 

F. THE CROSS CORRELATION FUNCTION 

The CCF has limited use for proving a signal or waveform to be LPI, but its 
analysis gives a very definite indication of the amount of commonality occurring between 
two different signals. Consider two different LPI signals used by two or more different 
netted LPI radars. The orthogonal polyphase signals being studied have the same 
parameters with respect to the waveform, power level, center frequency, bandwidth, code 
period, but the one difference is their phase code sequence. It is necessary for the LPI 
radars to maintain integrity of their signals in order to properly analyze the return signals 
for possible target detections. A zero cross-correlation shows that one radar’s signal will 
not interfere with the other. Using the equation below, a transmitted signal M(t) and the 

complex conjugate of a second signal v(t) are multiplied together and integrated over the 
length of the sequences: 

Tr 

Kv{'^)= \ + (15) 
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For CW polyphase signals, the integration is over the length of the code period. It is 
desirable for the result of the cross-correlation of the two polyphase signals to be zero. 
This represents an absence of interaction between the two signals. This will be illustrated 
in the following chapters of this paper [10]. 

G. SUMMARY 

In summary, the basic principles of signal analysis have been explained. In the 
next chapter, these principles will be applied to illustrate the performance of the OPS, 
Optimized OPS, and discrete frequency signals. In addition, the cross-correlation 
function will be used to show the non-interference characteristics of the orthogonal 
signals. 
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III. LPI WAVEFORMS FOR NETTED RADAR 


A. ORTHOGONAL POLYPHASE CODES 

Polyphase waveforms have long been employed in LPI radar systems due to the 
difficulty they introduce in adversarial ESR detection and analysis. Polyphase waveform 
sequences are identified as a series of phase sequences of finite length (), equal sub¬ 
code width (sub-code period, 4 ). For each phase sequence, the phase of each sub-code 

varies in a specific pattern over the length of the sequence. The transmitted signal and 
antenna bandwidth are not increased by an increase in sequence length; bandwidth (B) is 
directly related to the carrier frequency and inversely related to the number of cycles 

per phase (cpp) as shown by Equation(16). The bandwidth is also the inverse of the sub¬ 
code period. Code sequences can be designed to be longer to increase their processing 
gain and their resistance to compromise by adversary ESR without affecting the 
bandwidth size, although a longer code sequence will necessitate a more complex 
matched filter [ 1 ]. 

B = ^ = - (16) 

cpp tb 

We assume that the orthogonal polyphase code set consists of L signals with each 
signal containing Nc sub-codes as: 

{n) = = 1,2,...A^,/ = 1,...L 

The orthogonal polyphase waveforms are generally identified as a set (S) with three 
unique variables. The first variable is the code length ; this length is variable from 

two to greater than 1000. As the sequence length increases, so does radar performance in 
terms of auto-correlation side lobe levels decreasing. This performance increase with 
length is similar to that of other polyphase waveforms to include Frank codes and the PI - 
P4 codes [3]. As Nc increases, so does the complexity of the receiver, but the cost of 
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receiver design may outweigh the cost of adversary ESR improvements required to 
intercept and analyze the longer waveform. Computation time is increased drastically 
with values of Nc greater than 50 [3]. 

The second variable is the set size (L), usually from two to eight, but can be larger 
if required by design. The set size L is the number of codes of length Nc in a set. The 
value of L as applied to netted LPI radars determines the number of radar sites in the 
multisite radar system. The main complication with an increase value of L is the complex 
computation required to produce the final orthogonal polyphase waveform set S. For the 
same code length Nc it has been shown that the radar performance does not degrade as the 
set size L is increased [3]. 


The third variable to describe an orthogonal polyphase waveform is the number of 
admissible phase values (M) in the code sequence. The values are determined by the 
following general equation: 

,, . 271 r^2K 2 ; t , , , 

(j){n) e {0, —,2 —-1)—} = (17) 

MM M 

Using this equation for a value of 4, the choice of phase values would be: 
TT StT 

Q,—,71,and— and as the value of admissible phase values M increases so does the 


complexity of the receiver. Equation (18) illustrates a polycode phase set in an L by Nc 
phase matrix [3]. 


S{L,N^,M) 




(18) 


A(1)’^l(2),<2^l(3),..4(^c) 


Once the values of L, Nc, and M are established the next step in the process is to 
determine the code sets in the form of a phase matrix. The matrix will consist of L sets of 
phase sequences of length Nc- The method to determine the each phase value in the phase 
sequence is called simulated annealing (SA) algorithm. SA is required in the 
determination of multiple long phase sequences with low cross-correlation because 
algebraic determination of greater than two sequences is very difficult. SA was initially 
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designed for VLSI circuit design. The algorithm is an optimization algorithm with the 
benefit of not becoming fixed on local optima during the search process [3]. 


The optimization process is based on the comparison of two properties of the 
OPS, the autocorrelation and the cross correlation. Specifically the code set must satisfy 
or nearly satisfy the ACT and CCF shown in discrete form in Equations (19) and(20). 


A(s„k) = 


N-k 




^ Si(n)s*(n + k) = 0,0 < k < 


n=l 

N 


— ^ Si(n)s*(n + k) = 0,-N^ <k <0 


l = l,2,...L 


(19) 


n=-k+\ 


C{s^,s^,k) 


j^"Y^s^{n)s;{n + k) = 0,0 <k<N^ 

= ^ "''I p^q,l = \,2,...L 

Z Sp(n)sJ(n + k) = 0,-N^<k<0 

n=-k+l 


( 20 ) 


The basic cost function is given by Equation(21). The variable T is a weighting 
coefficient between the ACE and the CCF in the cost function. The principle of the 
function is to minimize the ACF sidelobe peaks and the CCF peaks. During the process 
using Equation (21) it is possible to obtain larger autocorrelation sidelobe levels with 
increased code length. Because of this a better optimization technique is to minimize 
total ACF sidelobe energy and CCF energy. This is shown in Equation(22). 

Ej =^max|A(<zi,,k)| + ;i ^ max\C(^^,^^,k)\,k ^0 (21) 

l=l p=l q=p+\ 


Z Z \C(t4,.k)\ (22) 

;=1 k=\ p=\ q=p+\k=-(N,-\) 

Given the predetermined values of L, M, and Nc, the minimization of Equation (22) using 
SA will generate a group of polyphase sequences that are automatically constrained by 
Equation (19) and Equation(20) [3]. 

The four basic steps of the SA process are as follows: 

• Step 1: Choose initial values of L, M, and Nc. 
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• Step 2: Choose annealing temperature. The initial starting temperature is 
based on the standard deviation of the initial cost distribution. This 
distribution is determined by sampling the number set produced by Equation 
(17) and illustrated by Equation( 18). 

• Step 3: Start the annealing process by “perturbing” the phase coding by 
randomly choosing a coding sequence from the waveform set. If the cost is 
reduced by the perturbation it will be accepted, if it increases it maybe 
accepted if it is within a probability tolerance based on the cost change and the 
current annealing temperature. This process is repeated until the standard 
deviation of the cost values becomes stable at the annealing temperature. 
Then the annealing temperature is changed until the standard deviation of the 
cost values becomes stable again. 

• Step 4: Stop the process when the cost is not reduced over the period of three 
consecutive temperature reductions [3]. 

1. Orthogonal Polyphase Sequence of 40 Suh-codes 
a. Waveform Description 

The orthogonal polyphase sequence code with 40 sub-codes (OPS 40) was 
designed with the following initial values, Nc = 40, L = 4, and M = 4. There are L = 4 
code sets associated with this particular design. The M = 4 available phases 

areO,——. The MATLAB® [21 code allows for user selection of all variables in 
2 2 

the Table 1. To prevent confusion and aid the user in creation analysis of signals a 
naming convention was established for the OPS 40 code. The OPS 40 code file names 
begin with “O 40 4 4” (OPS “N” “L” “M”). They are followed by the code set number 
(i.e. Cl), the carrier frequency (in kHz), the sampling frequency (in kHz), the code 
length, the number of cycles per phase, and finally the SNR (in dB, or “s” for signal 
only). 
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The OPS 40 code is selected as it is short enough that the algorithm used 
to generate it does not require a significant computation time, yet it is long enough that its 
performance exceeds that of shorter code lengths. The OPS 40 code set is displayed in 
Table 2. Each code sequence is listed vertically (L = 4) and the sub-code index is listed 
horizontally (Nc = 40) to show the difference in phase value between each sub-code in 
each sequence. The multiplier value to get the phase code value in radians is also shown 
in Table 2. This increased diversity provides an inherent invulnerability to the non- 
cooperative ESR and due to the nature of polyphase codes does not increase the signal 
bandwidth that is governed by the bit sub-code period [7]. The OPS 40 codes were 
generated for testing with the following characteristics shown in Table 1: 


Specific Signal Identification: 0 40 4 4 C1-1-7-40-1-S 

Carrier frequency ( /^) 

1000 Hz 

Sampling frequency (/^) 

7000 Hz 

Number of phases 

40 

SNR 

Signal Only 

Number of cycles per phase 

1 

Number of potential signals in the set 

4 


Table 1. OPS 40 Signal Characteristics. 
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Code 1 

Code 2 

Code 3 

Code 4 

Index 

XXTZ 

Rad 

XXTT 

Rad 

XXTT 

Rad 

XXTT 

Rad 

1 

0.50 

1.57 

1.50 

4.71 

1.50 

4.71 

0.50 

1.57 

2 

0.00 

0.00 

1.00 

3.14 

0.00 

0.00 

0.50 

1.57 

3 

0.00 

0.00 

1.00 

3.14 

0.50 

1.57 

0.00 

0.00 

4 

1.50 

4.71 

1.00 

3.14 

0.50 

1.57 

0.00 

0.00 

5 

0.50 

1.57 

0.50 

1.57 

1.50 

4.71 

0.00 

0.00 

6 

1.50 

4.71 

0.00 

0.00 

0.50 

1.57 

0.50 

1.57 

7 

1.50 

4.71 

0.00 

0.00 

1.50 

4.71 

1.50 

4.71 

8 

0.50 

1.57 

0.50 

1.57 

1.00 

3.14 

1.00 

3.14 

9 

0.50 

1.57 

0.50 

1.57 

0.50 

1.57 

1.50 

4.71 

10 

1.00 

3.14 

1.50 

4.71 

0.50 

1.57 

1.00 

3.14 

11 

0.50 

1.57 

0.00 

0.00 

0.50 

1.57 

0.00 

0.00 

12 

0.50 

1.57 

0.50 

1.57 

0.50 

1.57 

1.50 

4.71 

13 

1.50 

4.71 

1.00 

3.14 

0.50 

1.57 

1.50 

4.71 

14 

1.50 

4.71 

0.50 

1.57 

1.00 

3.14 

0.00 

0.00 

15 

1.00 

3.14 

1.00 

3.14 

1.00 

3.14 

1.00 

3.14 

16 

0.50 

1.57 

1.00 

3.14 

0.50 

1.57 

0.00 

0.00 

17 

0.50 

1.57 

0.50 

1.57 

0.00 

0.00 

0.00 

0.00 

18 

1.50 

4.71 

1.50 

4.71 

0.50 

1.57 

0.50 

1.57 

19 

0.00 

0.00 

0.00 

0.00 

1.50 

4.71 

1.50 

4.71 

20 

1.00 

3.14 

0.50 

1.57 

1.00 

3.14 

1.50 

4.71 

21 

0.00 

0.00 

1.50 

4.71 

0.50 

1.57 

0.50 

1.57 

22 

1.50 

4.71 

0.50 

1.57 

1.50 

4.71 

1.50 

4.71 

23 

0.50 

1.57 

0.50 

1.57 

0.00 

0.00 

0.00 

0.00 

24 

0.50 

1.57 

1.50 

4.71 

1.00 

3.14 

0.50 

1.57 

25 

1.00 

3.14 

0.50 

1.57 

0.50 

1.57 

0.00 

0.00 

26 

1.50 

4.71 

0.00 

0.00 

0.00 

0.00 

1.50 

4.71 

27 

1.50 

4.71 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

28 

1.50 

4.71 

1.50 

4.71 

1.00 

3.14 

0.00 

0.00 

29 

0.50 

1.57 

0.50 

1.57 

0.00 

0.00 

0.50 

1.57 

30 

1.50 

4.71 

0.50 

1.57 

1.00 

3.14 

0.50 

1.57 

31 

1.50 

4.71 

1.00 

3.14 

1.00 

3.14 

1.00 

3.14 

32 

1.50 

4.71 

1.00 

3.14 

0.50 

1.57 

1.50 

4.71 

33 

1.00 

3.14 

1.50 

4.71 

1.00 

3.14 

0.50 

1.57 

34 

0.50 

1.57 

0.50 

1.57 

1.00 

3.14 

0.50 

1.57 

35 

1.00 

3.14 

1.50 

4.71 

1.50 

4.71 

0.00 

0.00 

36 

1.00 

3.14 

0.50 

1.57 

0.50 

1.57 

1.50 

4.71 

37 

1.00 

3.14 

1.50 

4.71 

0.00 

0.00 

0.50 

1.57 

38 

1.00 

3.14 

0.50 

1.57 

0.00 

0.00 

1.00 

3.14 

39 

1.50 

4.71 

1.00 

3.14 

1.50 

4.71 

1.00 

3.14 

40 

0.50 

1.57 

0.50 

1.57 

0.00 

0.00 

1.00 

3.14 


Table 2. OPS Sequences with A^c = 40, L = 4, M = 4 (from [3]). 
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b. Waveform Characteristics from Simulation 


Figure 1 shows the power spectral density (PSD) versus frequency of the 
OPS 40 modulation code set #1 with Nc = 40. Figure 2 shows the discrete phase values 
for the OPS 40 modulation code set #1 with Nc = 40. The four phase values are readily 
apparent. Figure 3 shows the complex signal phase shift for the OPS 40 modulation code 
set #1. This plot best illustrates the different signal characteristics between the different 
code sets. The carrier frequency is /^ = 1 kHz, /^ = 7 kHz, and cpp = 1. All three plots 
were generated with the LPIT. 
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Figure 1. Power Spectral Density of OPS Sequence 

with Nc = 40, L = 4, M = 4 Code Set 1 (Signal only). 
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i - index for phase change 


Figure 2. OPS Sequence with = 40, L = 4, M = 4 Code Set 1 
showing one period of phase shifts (radians). 
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Figure 3. OPS Sequence with Nc = 40, L = 4, M = 4 Code Set 1 
showing five periods of phase shifts (degrees). 
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Figure 4 shows the PSD versus frequency of the OPS 40 modulation code 
set #2 with Nc = 40. Figure 5 shows the discrete phase values of the OPS 40 modulation 
code set #2 with Nc = 40. The four phase values are readily apparent. Figure 6 shows the 
“I + jQ” phase shift for the OPS 40 modulation code set #2. This plot best illustrates the 
different signal characteristics between the different code sets. The carrier frequency is 
/^ = 1 kHz, /^ = 7 kHz, and cpp = 1. All three plots were generated with the LPIT. 


10 

5 

0 


"D 

1 -10 

O) 

OJ 

E -15 

2 

S -20 

Q. 

CO 

I -25 
o 
CL 

-30 
-35 
-40 

0 500 1000 1500 2000 2500 3000 3500 

Frequency 

Figure 4. Power Speetral Density of OPS Sequenee 

with Nc = 40, L = 4, M = 4 Code Set 2 (Signal only). 
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Figure 5. OPS Sequence with A^c = 40, L = 4, M = 4 Code Set 2 
showing one period of phase shifts (radians). 
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Figure 6. OPS Sequence with Nc = 40, L = 4, M = 4 Code Set 2 
showing five periods of phase shifts (degrees). 

22 





































































Figure 7 shows the PSD versus frequency of the OPS 40 modulation code 
set #3 with Nc = 40. Figure 8 shows the discrete phase values for the OPS 40 modulation 
code set #3 with Nc = 40. The four phase values are readily apparent. Figure 9 shows the 
“I + jQ” phase shift for the OPS 40 modulation code set #3, this plot best illustrates the 
different signal characteristics between the different code sets. The carrier frequency is 
/^ = 1 kHz, /^ = 7 kHz, and cpp = 1. All three plots were generated with the LPIT. 
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Figure 7. Power Spectral Density of OPS Sequence 

with Nc = 40, L = 4, M = 4 Code Set 3 (Signal only). 
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Figure 8. OPS Sequence with A^c = 40, L = 4, M = 4 Code Set 3 
showing one period of phase shifts (radians). 
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Figure 9. OPS Sequence with Nc = 40, L = 4, M = 4 Code Set 3 
showing five periods of phase shifts (degrees). 
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Figure 10 shows the PSD versus frequency of the OPS 40 modulation 
code set #4 with Nc = 40. Figure 11 shows the discrete phase values for the OPS 40 
modulation code set #4 with Nc = 40. The four phase values are readily apparent. Figure 
12 shows the “I + jQ” phase shift for the OPS 40 modulation code set #4, this plot best 
illustrates the different signal characteristics between the different code sets. The carrier 
frequency is /^ = 1 kHz, /^ = 7 kHz, and cpp = 1. All three plots were generated with the 
LPIT. 
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Figure 10. Power Spectral Density of OPS Sequence 

with Nc = 40, L = 4, M = 4 Code Set 4 (Signal only). 


25 

























Figure 11. OPS Sequence with A^c = 40, L = 4, M = 4 Code Set 4 
showing one period of phase shifts (radians). 
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Figure 12. OPS Sequence with = 40, L = 4, M = 4 Code Set 4 
showing five periods of phase shifts (degrees). 
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Figure 13 (top plot) shows the normalized ACF and the corresponding 
sidelobe structure for the OPS 40 code set #1. The offset t is also normalized by the 
sub-code period 4 . The peak sidelobe level (PSL) can be seen at approximately -16 dB 
from the peak [3], Figure 13 (bottom plot) shows the PACF. Figure 14 shows the PAF. 
Based on the signal parameters = 1 kHz, /^ = 7 kHz, and cpp = 1 and using Equation 


(23) the PAF repeats every code period as shown by [1]: 


= A 


cpp* — 
V /'=' 


= 280 samples 


(23) 


All three plots were generated with ambfnV [16]. 



Figure 13. OPS Sequence with Ac = 40, L = 4, M = 4 Code Set 1 ACF (top) 

and PACF (bottom). 
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Figure 14. OPS Sequence with A^c = 40, L = 4, M = 4 Code Set 1 

Periodic Ambiguity Function Plot. 

To illustrate the non-interference nature of the OPS 40 signal the CCF was 
calculated for all six possible combinations of signals that would be utilized by a 
networked LPI radar system with four nodes. Each would be transmitting its own code 
set. Theoretically there would be no interference between each different code sets using 
similar signal parameters. Figure 15 shows the CCF for the OPS 40 code set #1 and code 
set #2 on the left contrasted with the ACF of OPS 40 code set #1 on the right. Figure 16 
shows the CCF for the OPS 40 code set #2 and code set #3 on the left contrasted with the 
ACF of OPS 40 code set #2 on the right. Figure 17 shows the CCF for the OPS 40 code 
set #3 and code set #4 on the left contrasted with the ACF of OPS 40 code set #3 on the 
right. Figure 18 shows the CCF for the OPS 40 code set #4 and code set #1 on the left 
contrasted with the ACF of OPS 40 code set #4 on the right. Figure 19 shows the CCF 
for the OPS 40 code set #4 and code set #2 on the left contrasted with the ACF of OPS 40 
code set #4 on the right. Figure 20 shows the CCF for the OPS 40 code set #3 and code 
set #1 on the left contrasted with the ACF of OPS 40 code set #3 on the right. 
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Figure 15. 


Figure 16. 
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OPS 40 Codes 1 and 2 CCF Plot (Left) and Code 1 ACF Plot (Right). 
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OPS 40 Codes 2 and 3 CCF Plot (Left) and Code 2 ACF Plot (Right). 
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Figure 17. 
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OPS 40 Codes 3 and 4 CCF Plot (Left) and Code 3 ACF Plot (Right). 


Figure 18. 
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OPS 40 Codes 4 and 1 CCF Plot (Left) and Code 4 ACF Plot (Right). 
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Figure 19. 


Figure 20. 
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OPS 40 Codes 4 and 2 CCF Plot (Left) and Code 4 ACF Plot (Right). 
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OPS 40 Codes 3 and 1 CCF Plot (Left) and Code 3 ACF Plot (Right). 
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2. Orthogonal Polyphase Sequence of 128 Suh-codes 


a. Waveform Description 

The orthogonal polyphase sequenee code with 128 sub-codes (OPS 128) 
was designed with the following initial values, Nc = 12S, L = 3, and M = 4. There are 
three code sets associated with this particular design. The four available phases 

areO,——. The MATLAB® 121 code allows for user selection of all variables in 
2 2 

the Table 3. To prevent confusion and aid the user in creation and analysis of signals a 
naming convention was established for the OPS 128 code. The entire OPS 128 code file 
names begin with O 128 3 4 (O “N” “L” “M”). They are followed by the code set 
number (i.e. Cl), the carrier frequency (in kHz), the sampling frequency (in kHz), the 
code length, the number of cycles per phase, and finally the SNR (in dB, or “s” for signal 
only). 

The OPS 128 code is much longer than the OPS 40 code. There is one 
less code set (L=3 vice L=4) reducing the computation time for signal design. It is much 
longer and its performance far exceeds that of shorter code lengths. The OPS 128 code is 
displayed in Table 4. Each code sequence is listed vertically (L = 3) and the sub-code 
index is listed horizontally (Nc = 128) to show the difference in phase value between each 
sub-code in each sequence. This increased diversity provides an increased invulnerability 
to non-cooperative ESR and due to the nature of polyphase codes does not increase the 
signal bandwidth that is governed by the sub-code period [7]. 
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The OPS 128 codes were generated for testing with the following 
characteristics shown in Table 3: 


Specific Signal Identification: 0 128 3 4 C1-1-7-128-1-S 

Carrier frequency (/^) 

1000 Hz 

Sampling frequency ( /^) 

7000 Hz 

Number of phases 

128 

SNR 

Signal Only 

Number of cycles per phase 

1 

Number of potential signals in the set 

3 


Table 3. OPS 128 Signal Characteristics. 
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Code 1 

Code 2 

Code 3 


Code 1 

Code 2 

Code 3 


Code 1 

Code 2 

Code 3 

11 

x*pi 

Rad 

x*pi 

Rad 

x*pi 

Rad 

B 

x*pi 

Rad 

x*pi 

Rad 

x*pi 

Rad 

# 

x*pi 

Rad 

x*pi 

Rad 

x*pi 

Rad 

1 

1.5 

4.71 

lilil 

EB 

1.5 

4.71 

EQ 

1.0 

3.14 

0.0 

0.00 

1.0 

3.14 

85 

0.0 

EB 

1.5 

4.71 

1.5 

4.71 

2 

1.5 

4.71 

0.5 

1.57 

1.0 

3.14 

B 

1.0 

3.14 

0.0 

0.00 

0.0 

EB 

86 

0.0 

EB 

0.5 

1.57 

1.5 

4.71 

H 

1.5 

4.71 

1.0 

3.14 

0.0 

EB 

B 

0.0 

0.00 

0.0 

0.00 

1.0 

3.14 

87 

0.5 

1.57 

1.0 

3.14 

0.0 

0.00 


1.5 

4.71 

1.0 

3.14 

0.0 

0.00 

B 

1.5 

4.71 

1.5 

4.71 

1.5 

4.71 

88 

0.5 

1.57 

0.5 

1.57 

0.5 

1.57 

B 

0.0 

EB 

0.0 

EB 

1.0 

3.14 

B 

1.5 

4.71 

1.0 

3.14 

0.0 

EB 

89 

1.0 

3.14 

0.5 

1.57 

1.0 

3.14 


0.0 

0.00 

1.5 

4.71 

1.5 

4.71 

B 

1.0 

3.14 

1.5 

4.71 

0.0 

0.00 

90 

0.5 

1.57 

1.0 

3.14 

0.0 

0.00 

H 

0.0 

0.00 

1.5 

4.71 

0.0 

0.00 

B 

0.0 

0.00 

1.0 

3.14 

1.0 

3.14 

91 

1.0 

3.14 

1.5 

4.71 

1.5 

4.71 

8 

1.5 

4.71 

0.5 

1.57 

0.0 

0.00 

B 

1.5 

4.71 

1.5 

4.71 

0.5 

1.57 

92 

0.5 

1.57 

0.0 

0.00 

1.5 

4.71 

9 

0.0 

0.00 

1.5 

4.71 

0.0 

0.00 

B 

1.0 

3.14 

1.5 

4.71 

1.0 

3.14 

93 

1.5 

4.71 

0.5 

1.57 

1.0 

3.14 

IQ 

1.5 

4.71 

0.5 

1.57 

0.0 

0.00 

m 

1.5 

4.71 

0.0 

0.00 

0.0 

0.00 

94 

0.0 

0.00 

0.0 

0.00 

1.5 

4.71 

in 

0.0 

EB 

1.5 

4.71 

1.0 

3.14 

B 

1.5 

4.71 

1.0 

3.14 

1.5 

4.71 

95 

0.5 

1.57 

0.5 

1.57 

1.5 

4.71 

m 

1.5 

4.71 

lilil 

EB 

1.0 

3.14 

m 

0.0 

0.00 

0.0 

0.00 

0.0 

EB 

96 

1.0 

3.14 

1.5 

4.71 

1.0 

3.14 

13 

0.0 

0.00 

0.5 

1.57 

1.5 

4.71 

B 

1.0 

3.14 

1.0 

3.14 

1.0 

3.14 

97 

0.5 

1.57 

lilil 

EB 

1.5 

4.71 

14 

0.0 

0.00 

0.0 

EB 

0.5 

1.57 

B 

1.0 

3.14 

0.5 

1.57 

0.0 

EB 

98 

lilil 

EB 

0.5 

1.57 

1.0 

3.14 

IQ 

1.5 

4.71 

liU 

1.57 

1.5 

4.71 

B 

0.0 

0.00 

0.5 

1.57 

0.5 

1.57 

99 

1.0 

3.14 

1.5 

4.71 

0.5 

1.57 

IQ 

1.5 

4.71 

liU 

1.57 

1.5 

4.71 

B 

1.0 

3.14 

0.5 

1.57 

0.5 

1.57 

100 

1.0 

3.14 

1.5 

4.71 

1.5 

4.71 

IQ 

0.5 

1.57 

0.0 

EB 

0.0 

EB 

B 

0.5 

1.57 

1.6 

4.87 

1.0 

3.14 

101 

1.0 

3.14 

0.5 

1.57 

0.0 

0.00 

IQ 

1.5 

4.71 

0.0 

0.00 

0.0 

0.00 

B 

0.5 

1.57 

0.0 

0.00 

0.0 

0.00 

102 

1.5 

4.71 

1.5 

4.71 

0.5 

1.57 

IQ 

1.0 

3.14 

1.0 

3.14 

1.0 

3.14 

B 

0.5 

1.57 

1.5 

4.71 

1.0 

3.14 

103 

0.0 

0.00 

1.5 

4.71 

0.5 

1.57 

Bil 

0.0 

0.00 

1.0 

3.14 

0.5 

1.57 

B 

1.5 

4.71 

1.0 

3.14 

1.5 

4.71 

104 

0.0 

0.00 

0.0 

0.00 

1.0 

3.14 

Bl 

0.0 

0.00 

1.5 

4.71 

0.5 

1.57 

B 

0.0 

0.00 

1.5 

4.71 

1.0 

3.14 

105 

1.5 

4.71 

1.0 

3.14 

0.5 

1.57 


0.0 

0.00 

1.5 

4.71 

0.0 

0.00 

B 

1.5 

4.71 

0.0 

0.00 

0.5 

1.57 

106 

1.5 

4.71 

1.0 

3.14 

0.5 

1.57 

Bel 

0.0 

EB 

1.0 

3.14 

0.5 

1.57 

B 

0.5 

1.57 

1.5 

4.71 

0.0 

EB 

107 

0.5 

1.57 

1.5 

4.71 

0.5 

1.57 

B 

1.5 

4.71 

0.0 

EB 

1.5 

4.71 

B 

1.0 

3.14 

1.0 

3.14 

1.5 

4.71 

108 

0.5 

1.57 

1.5 

4.71 

0.5 

1.57 

B 

0.5 

1.57 

1.5 

4.71 

0.0 

EB 

B 

0.0 

0.00 

1.0 

3.14 

1.5 

4.71 

109 

lilil 

EB 

0.5 

1.57 

0.5 

1.57 

B 

1.0 

3.14 

0.5 

1.57 

1.0 

3.14 

B 

1.0 

3.14 

0.5 

1.57 

1.0 

3.14 

110 

lilil 

1.57 

1.0 

3.14 

1.5 

4.71 

B 

1.5 

4.71 

0.0 

EB 

0.0 

EB 

B 

1.5 

4.71 

0.0 

0.00 

1.0 

3.14 

111 

0.0 

EB 

0.0 

EB 

1.0 

3.14 

B 

0.5 

1.57 

1.5 

4.71 

1.0 

3.14 

m 

1.0 

3.14 

0.0 

0.00 

0.5 

1.57 

112 

0.0 

0.00 

1.0 

3.14 

1.0 

3.14 

B 

0.5 

1.57 

1.0 

3.14 

1.0 

3.14 

B 

0.0 

0.00 

1.5 

4.71 

0.5 

1.57 

113 

1.0 

3.14 

1.0 

3.14 

1.5 

4.71 

B 

1.5 

4.71 

1.5 

4.71 

1.5 

4.71 

B 

1.0 

3.14 

1.0 

3.14 

0.0 

EB 

114 

1.0 

3.14 

1.5 

4.71 

0.5 

1.57 

B 

0.5 

1.57 

0.5 

1.57 

0.0 

0.00 

B 

1.0 

3.14 

0.0 

0.00 

1.0 

3.14 

115 

1.5 

4.71 

1.5 

4.71 

0.0 

0.00 

B 

1.5 

4.71 

0.0 

0.00 

0.0 

0.00 

B 

0.0 

0.00 

0.5 

1.57 

0.0 

0.00 

116 

1.5 

4.71 

1.5 

4.71 

1.0 

3.14 

B 

0.5 

1.57 

1.0 

3.14 

1.0 

3.14 

B 

0.5 

1.57 

1.0 

3.14 

0.0 

0.00 

117 

0.5 

1.57 

1.5 

4.71 

0.5 

1.57 

B 

0.0 

0.00 

1.0 

3.14 

0.5 

1.57 

B 

1.0 

3.14 

0.0 

0.00 

0.0 

0.00 

118 

0.5 

1.57 

0.0 

0.00 

1.5 

4.71 

B 

0.5 

1.57 

lilil 

EB 

1.0 

3.14 

B 

1.0 

3.14 

0.5 

1.57 

lilil 

EB 

119 

1.5 

4.71 

0.5 

1.57 

1.0 

3.14 

B 

0.5 

1.57 

lilil 

EB 

0.0 

EB 

B 

1.5 

4.71 

1.5 

4.71 

1.0 

3.14 

120 

0.0 

EB 

1.5 

4.71 

1.0 

3.14 

B 

1.0 

3.14 

lilil 

EB 

1.5 

4.71 

B 

1.0 

3.14 

1.0 

3.14 

1.5 

4.71 

121 

1.0 

3.14 

1.0 

3.14 

0.0 

0.00 

B 

lilil 

EB 

1.0 

3.14 

0.0 

EB 

B 

0.5 

1.57 

1.5 

4.71 

1.0 

3.14 

122 

0.0 

EB 

1.5 

4.71 

0.0 

0.00 

39 

0.0 

0.00 

1.0 

3.14 

1.0 

3.14 

B 

0.0 

0.00 

1.0 

3.14 

1.0 

3.14 

123 

0.5 

1.57 

1.5 

4.71 

1.0 

3.14 

B 

0.0 

0.00 

0.5 

1.57 

0.0 

0.00 

B 

0.5 

1.57 

0.0 

0.00 

0.5 

1.57 

124 

1.5 

4.71 

1.5 

4.71 

1.0 

3.14 

B 

0.5 

1.57 

0.5 

1.57 

0.5 

1.57 

B 

1.5 

4.71 

1.5 

4.71 

1.5 

4.71 

125 

0.0 

0.00 

1.5 

4.71 

1.5 

4.71 

B 

0.0 

0.00 

1.0 

3.14 

0.5 

1.57 

B 

0.0 

0.00 

0.5 

1.57 

1.0 

3.14 

126 

1.0 

3.14 

1.5 

4.71 

1.5 

4.71 


127 

0.5 

1.57 

1.0 

3.14 

1.5 

4.71 

128 

0.0 

0.00 

1.5 

4.71 

0.5 

1.57 


Table 4. OPS Sequences with Nc= 128, L = 3, M = 4 (from [3]). 
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b. Waveform Characteristics from Simulation 


Figure 21 shows the PSD versus frequency of the OPS 128 modulation 
code set #1 with Nc = 128. Figure 22 shows the discrete phase values for the OPS 128 
modulation code set #1 with Nc = 128. The four phase values are readily apparent. 
Figure 23 shows the “1 + jQ” phase shift for the OPS 128 modulation code set #1. This 
plot best illustrates the different signal characteristics between the different code sets. 
The carrier frequency is /^ = 1 kHz, /^ = 7 kHz, and cpp = 1. All three plots were 
generated with the LPIT. 


PSD of I Phase Shift & no Noise 



Figure 21. Power Spectral Density of OPS 128 Code 1 (Signal only). 
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Figure 22. OPS 128 Code 1 showing one period of phase shifts (radians). 



Figure 23. OPS 128 Code 1 showing five periods of phase shifts (degrees). 
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Figure 24 shows the PSD versus frequency of the OPS 128 modulation 
code set #2 with Nc = 128. Figure 25 shows the discrete phase values for the OPS 128 
modulation code set #2 with Nc = 128, the four phase values are readily apparent. Figure 
26 shows the “I + jQ” phase shift for the OPS 128 modulation code set #2. This plot best 
illustrates the different signal characteristics between the different code sets. The carrier 
frequency is /^ = 1 kHz, /^ = 7 kHz, and cpp = 1. All three plots were generated with the 
LPIT. 
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Figure 24. Power Spectral Density of OPS 128 Code 2 (Signal only). 
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Figure 25. OPS 128 Code 2 showing one period of phase shifts (radians). 
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Figure 26. OPS 128 Code 2 showing five periods of phase shifts (degrees). 
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Figure 27 shows the PSD versus frequency of the OPS 128 modulation 
code set #3 with Nc = 128. Figure 28 shows the discrete phase values for the OPS 128 
modulation code set #3 with Nc = 128, the four phase values are readily apparent. Figure 
29 shows the “I + jQ” phase shift for the OPS 128 modulation code set #3. This plot best 
illustrates the different signal characteristics between the different code sets. The carrier 
frequency is /^ = 1 kHz, /^ = 7 kHz, and cpp = 1. All three plots were generated with the 
LPIT. 
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Figure 27. Power Spectral Density of OPS 128 Code 3 (Signal only). 
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Figure 28. OPS 128 Code 3 showing one period of phase shifts (radians). 



Figure 29. OPS 128 Code 3 showing five periods of phase shifts (degrees). 
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Figure 30 (top plot) shows the ACF and the corresponding sidelobe 
structure for the OPS 128 code set #1. The peak sidelobe level (PSL) can be seen at 
approximately -21 dB from the peak [3]. Figure 30 (bottom plot) shows the PACF. 
Figure 31 shows the PAF. Based on the signal parameters /^ = 1 kHz, /^ = 7 kHz, and 


cpp = 1 and using Equation (24) the PAF repeats every code period as shown by [1]: 

= 896 samples (24) 


^ fs^ 

cpp*^ 


b =N 

SC c ^ r 

V Ay 

All three plots were generated with ambfin? [16], 




Figure 30. OPS 128 Code 1 ACF (top) and PACF (bottom). 
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Figure 31. OPS 128 Code 1 Periodie Ambiguity Function Plot. 

To illustrate the non-interference nature of the OPS 128 signal the CCF 
was calculated for all three possible combinations of signals that would be utilized by a 
networked LPl radar system with three nodes. Each would be transmitting its own code 
set. Theoretically there would be no interference between each different code sets using 
similar signal parameters. Figure 32 shows the CCF for the OPS 128 code set #1 and 
code set #2 on the left contrasted with the ACF of OPS 40 code set #1 on the right. 
Figure 33 shows the CCF for the OPS 128 code set #2 and code set #3 on the left 
contrasted with the ACF of OPS 40 code set #2 on the right. Figure 34 shows the CCF 
for the OPS 128 code set #3 and code set #1 on the left contrasted with the ACF of OPS 
40 code set #3 on the right. 
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Figure 32. 


OPS 128 Codes 1 and 2 CCF Plot (Left) and Code 1 ACF Plot (Right). 


Figure 33. 
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OPS 128 Codes 2 and 3CCF Plot (Left) and Code 2 ACF Plot (Right). 
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Figure 34. OPS 128 Codes 3 and 1 CCF Plot (Left) and Code 3 ACF Plot (Right). 

3. Optimized Orthogonal Polyphase Sequences with Doppler Tolerance 

The Optimized OPS by Kahn is an improvement on Deng’s orthogonal polyphase 
sequenees and are obtained through the use of an improved design tool known as eross 
entropy (CE). Using this optimization teehnique, the new polyphase sequenees provide 
an improved Doppler toleranee in addition to being orthogonal. Deng’s OPS required 
Doppler eorreetion by the radar reeeiver while the Khan eode is designed to be inherently 
Doppler resistant [5]. 


cross-corr with 
O 128 3 4 CM-7-128-1-s 
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Offset 


ACF 

O 128 3 4C3-1-7-128-1-S 



0 

Offset 


a. Waveform Description 

The Optimized OPS 40 code was designed with the following initial 

values, L = 40, M = 3, and y = 4. There are three code sets associated with this particular 

design. There are many more available phases for use in the Optimized OPS. These 

phases are more complex than for OPS 40, but this complexity provides the Doppler 

tolerance characteristic of the signal design. The selection of code phases is determined 
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from an array based on the value of Y. The MATLAB [2] code allows for user selection 
of all variables in Table 5. To prevent confusion and aid the user in creation of signals a 
naming convention was established for the Optimized OPS 40 code. The entire 
Optimized OPS 40 code fde names begin with KO 40 4 4 (KO “N” “M” “y”). They are 
followed by the set number (i.e. Cl), the carrier frequency (in kHz), the sampling 
frequency (in kHz), the code length, the number of cycles per phase, and finally the SNR 
(in dB, or “s” for signal only). 

The Optimized OPS 40 code is selected as the length is short enough that 
the algorithm does not require a significant computation time, yet it is long enough that 
its performance exceeds that of shorter code lengths. The Optimized OPS 40 code is 
displayed in Table 6. Each code sequence is listed vertically (M = 3) and the sub-code 
index is listed horizontally (Nc = 40) to show the difference in phase value between each 
sub-code in each sequence. Equation (25) shows the method of determining the 
polyphase code value, where m is the number of code sequences (or radar transmitters), P 
is the number of phases to choose from, and y is the number of polyphase code value sets 
to select from on either side of the center of possible sets illustrated below: 

=|(|,-,-| + l)^(0,l,2,...,(P-2))| (25) 

The center-most (i=4) phase value sets to choose from using Equation (25) and y=4 and 
P=10 are from the following set{2;r/40(0,1,2,...,8)}. Figure 35 illustrates the tendency 

for values of the polyphase codes available to get larger as the value of i moves away 
from the center value. 

This increased diversity in phase values provides an inh erent 
invulnerability to non-cooperative ESR and due to the nature of polyphase codes does not 
increase the signal bandwidth that is governed by the sub-code period. The harmonic 
relationship between polyphase sequences from one phase value to the next is the 
property that provides the Doppler tolerance. Building this harmonic nature into the 
creation of the polyphase code sequence decreases the loss due to the effect of a Doppler 
shift [7]. The Optimized OPS 40 code was generated for testing with the following 
characteristics as shown in Table 5. 
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Specific Signal Identification: KO 40 3 4 C1-1-7-40-1-s 

Carrier frequency ( /^) 

1000 Hz 

Sampling frequency (/^) 

7000 Hz 

Number of phases 

40 

SNR 

Signal Only 

Number of cycles per phase 

1 

Number of potential signals in the set 

3 


Table 5. Optimized OPS 40 Signal Characteristics. 
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Code 1 

Code 2 

Code 3 

Index 

XXTl 

Rad 

XXTl 

Rad 

XXTl 

Rad 

1 

1.60 

5.03 

0.60 

1.88 

0.80 

2.51 

2 

1.00 

3.14 

0.20 

0.63 

0.40 

1.26 

3 

1.20 

3.77 

1.00 

3.14 

0.00 

0.00 

4 

0.80 

2.51 

1.20 

3.77 

0.00 

0.00 

5 

1.80 

5.65 

1.80 

5.65 

1.00 

3.14 

6 

0.40 

1.26 

1.80 

5.65 

1.80 

5.65 

7 

1.35 

4.24 

0.15 

0.47 

1.05 

3.30 

8 

1.05 

3.30 

0.60 

1.88 

0.45 

1.41 

9 

0.30 

0.94 

1.35 

4.24 

1.20 

3.77 

10 

0.15 

0.47 

0.30 

0.94 

0.45 

1.41 

11 

1.35 

4.24 

1.35 

4.24 

1.05 

3.30 

12 

1.35 

4.24 

1.20 

3.77 

1.35 

4.24 

13 

0.20 

0.63 

0.50 

1.57 

0.80 

2.51 

14 

0.00 

0.00 

0.80 

2.51 

0.90 

2.83 

15 

0.10 

0.31 

0.50 

1.57 

0.10 

0.31 

16 

0.80 

2.51 

0.00 

0.00 

0.20 

0.63 

17 

0.00 

0.00 

0.70 

2.20 

0.90 

2.83 

18 

0.60 

1.88 

0.80 

2.51 

0.90 

2.83 

19 

0.05 

0.16 

0.00 

0.00 

0.00 

0.00 

20 

0.20 

0.63 

0.05 

0.16 

0.05 

0.16 

21 

0.05 

0.16 

0.10 

0.31 

0.05 

0.16 

22 

0.25 

0.79 

0.05 

0.16 

0.25 

0.79 

23 

0.90 

2.83 

0.80 

2.51 

0.00 

0.00 

24 

0.90 

2.83 

0.10 

0.31 

0.90 

2.83 

25 

0.20 

0.63 

0.80 

2.51 

0.00 

0.00 

26 

0.70 

2.20 

0.10 

0.31 

0.80 

2.51 

27 

0.00 

0.00 

0.90 

2.83 

0.80 

2.51 

28 

0.20 

0.63 

0.60 

1.88 

0.20 

0.63 

29 

1.20 

3.77 

1.20 

3.77 

1.35 

4.24 

30 

1.35 

4.24 

0.45 

1.41 

0.15 

0.47 

31 

0.00 

0.00 

0.45 

1.41 

1.35 

4.24 

32 

0.00 

0.00 

0.00 

0.00 

1.35 

4.24 

33 

1.20 

3.77 

1.35 

4.24 

0.30 

0.94 

34 

0.30 

0.94 

1.20 

3.77 

0.75 

2.36 

35 

1.40 

4.40 

0.40 

1.26 

1.20 

3.77 

36 

0.60 

1.88 

1.20 

3.77 

0.60 

1.88 

37 

0.20 

0.63 

1.40 

4.40 

0.60 

1.88 

38 

1.20 

3.77 

0.60 

1.88 

1.20 

3.77 

39 

1.20 

3.77 

0.00 

0.00 

1.40 

4.40 

40 

1.60 

5.03 

1.80 

5.65 

1.40 

4.40 


Table 6. Optimized OPS Sequences L = 40, M = 3, Y=4 (from [5]). 
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Figure 35. Graphical Representation of Potential Phase Code Values. 

b. Waveform Characteristics from Simulation 

Figure 36 shows the PSD versus frequency of the Optimized OPS 40 
modulation code set #1 with Nc = 40. Figure 37 shows the discrete phase values for the 
Optimized OPS 40 modulation code set #1 with Nc = 40, the four phase values are readily 
apparent. Figure 38 shows the “I + jQ” phase shift for the Optimized OPS 40 modulation 
code set #1. This plot best illustrates the different signal characteristics between the 
different code sets. The carrier frequency is f^ =lkHz, =7 kHz, and cpp = 1. All 
three plots were generated with the LPIT. 
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Figure 36. 


Figure 37. 
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Power Spectral Density of Optimized OPS 40 Code 1 (Signal only). 



Optimized OPS 40 Code 1 showing one period of phase shifts (radians). 


49 

































































0 


Phase Shift from l+jQ 


B 

0) 


CO 

CD 

C/) 

C3 

CL 


-10 

-20 

-30 

-40 

-50 

-60 

-70 

-80 

-90 

-100 



III. 

1-: 

1-1-^^- 

- 

^ 4 ^ 

4 4 4 


4 4 4 

4 T 4 





- 

4 4 4 

4 4 4 


' ' 4 4 4 

4 4 4 

‘ A La^ 




1 

_ 1 

1 _ 1 _ 1 _ 1 _i 

I_i 

[_i_^^_ 


0 20 40 60 80 100 120 140 160 180 200 


i - index for phase change 


Figure 38. Optimized OPS 40 Code 1 showing five periods of phase 

shifts (degrees). 

Figure 39 shows the PSD versus frequeney of the Optimized OPS 40 
modulation eode set #2 with Nc = 40. Figure 40 shows the diserete phase values for the 
optimized OPS 40 modulation eode set #2 with Nc = 40, the four phase values are readily 
apparent. Figure 41 shows the “I + jQ” phase shift optimized for the OPS 40 modulation 
eode set #2. This plot best illustrates the different signal eharaeteristies between the 
different eode sets. The earrier frequeney is =lkHz, =7 kHz, and cpp = 1. All 
three plots were generated with the LPIT. 
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Figure 39. Power Speetral Density of Optimized OPS 40 Code 2 (Signal only). 



Figure 40. Optimized OPS 40 Code 2 showing one period of phase 

shifts (radians). 
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Figure 41. Optimized OPS 40 Code 2 showing five periods of phase 

shifts (degrees). 

Figure 42 shows the PSD versus frequency of the Optimized OPS 40 
modulation code set #3 with Nc = 40. Figure 43 shows the discrete phase values for the 
optimized OPS 40 modulation code set #3 with Nc = 40, the four phase values are readily 
apparent. Figure 44 shows the “I + jQ” phase shift optimized for the OPS 40 modulation 
code set #3. This plot best illustrates the different signal characteristics between the 
different code sets. The carrier frequency is =lkHz, =7 kHz, and cpp = 1. All 
three plots were generated with the LPIT. 
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Figure 42. Power Speetral Density of Optimized OPS 40 Code 3 (Signal only). 



Figure 43. Optimized OPS 40 Code 3 showing one period of phase 

shifts (radians). 


53 



































































B 

0) 


CO 

CD 

C/) 

CO 

x: 

CL 


5 

0 

-5 

-10 

-15 

-20 

-25 

-30 

-35 

0 20 40 60 80 100 120 140 160 180 200 

i - index for phase change 


Phase Shift from l+jQ 



1 

J 

1 1 

a L 

1-: 

1-!-^- 

- 

U 


- ■ -1- 

- h -1- 


-1-h- 

- 1- 


-1- 

-h-^- 


-1-h- 

- 1-- 


n 

T r 


^. T 

r 

— 

1_1_ 

_ 1 _ 1 _i 

I_i 

1 U' ; 

I_^^_ 

_^_ 


Figure 44. Optimized OPS 40 Code 3 showing five periods of phase 

shifts (degrees). 


Figure 45 (top plot) shows the ACF and the eorresponding sidelobe 
strueture for the optimized OPS 40 eode set #1. The PSL ean be seen at approximately 
-17 dB from the peak in Figure 45 (bottom plot). All three plots were generated with 
ambfn? [16]. Figure 46 shows the PAF. Based on the signal parameters =lkHz, 


f,=l kHz, and cpp = 1 and using Equation (26) the PAF repeats every eode period as 


shown by [1]: 


b.. = A. 


cpp* — 
V 


= 280 samples 


(26) 
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Figure 45. Optimized OPS 40 Code 1 ACF (top) and PACF (bottom). 



Figure 46. Optimized OPS 40 Code 1 Periodic Ambiguity Function Plot. 
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To illustrate the non-interference nature of the Optimized OPS 40 signal 
the CCF was calculated for all three possible combinations of signals that would be 
utilized by a networked LPI radar system with three nodes. Each would be transmitting 
its own code set. Theoretically there would be no interference between each different 
code sets using similar signal parameters. Figure 47 shows the CCF for the optimized 
OPS 40 code set #1 and code set #2 on the left contrasted with the ACF of optimized 
OPS 40 code set #1 on the right. Figure 48 shows the CCF for the optimized OPS 40 
code set #2 and code set #3 on the left contrasted with the ACF of optimized OPS 40 
code set #2 on the right. Figure 49 shows the CCF for the optimized OPS 40 code set #3 
and code set #1 on the left contrasted with the ACF of optimized OPS 40 code set #3 on 
the right. 
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Figure 47. Opt. OPS 40 Codes 1 and 2 CCF Plot (Left) and Code 1 

ACF Plot (Right). 
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Figure 48. Opt. OPS 40 Codes 2 and 3 CCF Plot (Left) and Code 2 

ACF Plot (Right). 
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Figure 49. Opt. OPS 40 Codes 3 and 1 CCF Plot (Left) and Code 3 

ACF Plot (Right). 
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B. DISCRETE FREQUENCY CODING 


1. DEC 32 Frequency Series 

The orthogonal discrete frequency coding (DFC) waveforms are identified as a set 
of L frequency hopping sequences. Each frequency hopping sequence is a CW waveform 
defined by the number of contiguous frequency codes of each duration tp. As length 

Np increases, so does radar performance in terms of the auto-correlation side lobe levels 
decreasing. Computation time deriving the codes is increased with increasing values 
of Np. Similar frequency hopping waveforms have been developed such as the Costas 

arrays, but for use in a single radar system. Algebraic calculations for two or more 
different Costas arrays do not guarantee a low cross correlation between the different 
waveforms [3]. 

The next step is the formation of a frequency code matrix, which consists of L sets 
of frequency codes of length A^. SA is required in the determination of multiple 

frequency coded sequences with low cross correlation because algebraic determination of 
greater than two sequences is very difficult [3]. SA was initially designed for VLSI 
circuit design [3]. The algorithm is an optimization algorithm with the benefit of not 
becoming fixed on local optima during the search process. The optimization process is 
based on the comparison of two properties of the DFC signal, the autocorrelation and the 
cross correlation. Specifically the code set must satisfy Equation (27) and Equation(28) 
[5]. 

A(s„T)=-l-is,(t)s;(t-T)dtrY‘,^ ■ ( 27 ) 

[^0, otherwise 

, r) = -^ J (t)sj {t-T)dtKO,p^q,p,q = 1,2,.. .L (28) 

The cost function is given by Equation(29). The principle of the cost function is to 
minimize the autocorrelation sidelobe peaks and the cross correlation peaks. The 
optimization technique is to minimize the total autocorrelation sidelobe energy and cross 
correlation energy. This is shown in Equation(29). 
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( 29 ) 


l=\ 2- P=1 ^ = /?+l 2- 

Given the predetermined values of L and the minimization of Equation (29) utilizing 

SA will generate a group of frequency coded sequences that are automatically 
constrained by Equation (27) and Equation(28) [5]. For signal generation the frequency 
multiplier value was multiplied by a base frequency to achieve the 32 frequency values 
for the signal. Most of the DFC signals that were generated had a base frequency of 1 
kHz, this yielding a frequency range from 1 kHz to 32 kHz for the signal to be analyzed. 

a. Waveform Description 

The Discrete Frequency Coding sequence was generated for testing with 
the following characteristics shown in Table 7. Table 8 shows the frequency multiplier 
values for DFC signal with 32 frequencies. 


Specific Signal Identification: DFC-32-1-7-1000 

Base frequency (/ ) 

1000 Hz 

Sampling frequency (f) 

7000 Hz 

Number of frequencies (Np) 

32 

SNR 

Signal Only 

Frequency duration (tp) 

0.001 s 

Number of potential signals in the set 

3 


Table 7. DFC Signal Characteristics. 
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Code 1 

Code 2 

Code 3 

Index 

Frequency Multiplier 

1 

32 

3 

32 

2 

29 

12 

14 

3 

2 

13 

19 

4 

31 

15 

21 

5 

4 

30 

12 

6 

9 

17 

8 

7 

11 

7 

30 

8 

13 

2 

28 

9 

8 

10 

9 

10 

26 

22 

23 

11 

12 

25 

1 

12 

16 

24 

2 

13 

24 

6 

22 

14 

14 

27 

15 

15 

28 

20 

10 

16 

23 

4 

18 

17 

27 

8 

6 

18 

18 

31 

26 

19 

6 

14 

27 

20 

22 

9 

11 

21 

7 

21 

20 

22 

3 

18 

17 

23 

30 

19 

31 

24 

16 

5 

16 

25 

17 

1 

13 

26 

20 

29 

24 

27 

10 

32 

5 

28 

21 

16 

4 

29 

19 

28 

7 

30 

1 

11 

3 

31 

25 

26 

29 

32 

5 

23 

25 


Table 8. DFC Frequency-Coding Sequences (based on [4]). 

b. Waveform Characteristics from Simulation 

Figure 50 shows the power spectrum magnitude versus frequency of the 
DFC modulation code set #1 with Np = 32. The base frequency is / =lkHz, 

/j =100 kHz, and = 0.001 s. The plot was generated with the LPIT. 
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Figure 50. Power Spectral Density of DFC Set 1 (Signal Only). 


Figure 51 shows the power spectrum magnitude versus frequency of the 
DFC modulation code set #2 with Np = 32. The base frequency is / =lkHz, 

/j =100 kHz, and tp = 0.001 s. The plot was generated with the LPIT. 
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Figure 51. Power Spectral Density of DFC Set 2 (Signal Only). 


Figure 52 shows the power spectrum magnitude versus frequency of the 
DFC modulation code set #3 with Np = 32. The base frequency is / =lkHz, 

/j =100 kHz, and tp = 0.001 s. The plot was generated with the LPIT. 
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Figure 52. Power Spectral Density of DFC Set 3 (Signal Only). 


Figure 53 (top plot) shows the ACF and the corresponding sidelobe 
structure for the OPS 40 code set #1. The peak sidelobe level (PSL) is slightly higher 
than -20 dB from the peak [4]. Figure 53 (bottom plot) shows the PACF. Figure 54 
shows the PAF. The carrier frequency is / = IkHz, = 100 kHz, and tp = 0.001 s. All 
three plots were generated with ambfn? [16]. 


63 




















Figure 53. DFC 32 Code 1 ACF (top) and PACF (bottom). 



Figure 54. DFC 32 Code 1 Periodic Ambiguity Function Plot. 
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To illustrate the non-interfering nature of the DFC 32 signal the CCF was 
applied to all three possible combinations of signals that would be utilized by a 
networked LPI radar system with three nodes. Each radar would be transmitting its own 
code set. Theoretically there would be no interference between each different code sets 
using similar signal parameters. Figure 55 shows the CCF for the DFC 32 code set #1 
and code set #2 on the left contrasted with the ACF of DFC 32 code set #1 on the right. 
Figure 56 shows the CCF for the DFC 32 code set #2 and code set #3 on the left 
contrasted with the ACF of DFC 32 code set #2 on the right. Figure 57 shows the CCF 
for the DFC 32 code set #3 and code set #1 on the left contrasted with the ACF of DFC 
32 code set #3 on the right. 
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Figure 55. DFC 32 Codes 1 and 2 CCF Plot (Left) and Code 1 ACF Plot (Right). 
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Figure 56. 


Figure 57. 
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DFC 32 Codes 2 and 3 CCF Plot (Left) and Code 2 ACF Plot (Right). 
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DFC 32 Codes 3 and 1 CCF Plot (Left) and Code 3 ACF Plot (Right). 
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C. SUMMARY 

In summary, Table 9 shows the PSL for each of the codes examined in this 
chapter. In the next chapter, the codes generated will be analyzed using time-frequency 
and bi-frequency techniques to extract the signal parameters. 


Code 

N/Nf 

PSL (dB) 

OPS 

40 

-16 

OPS 

128 

-21 

Optimized OPS 

40 

-17 

DFC 

32 

-20 


Table 9. PSL of Orthogonal Codes. 


67 




THIS PAGE INTENTIONALLY LEFT BLANK 


68 



IV. WAVEFORM PARAMETER EXTRACTION TECHNIQUES 


A. WIGNER-VILLE DISTRIBUTION 

1. Introduction 

The first waveform analysis performed on the signals was the Wigner-Ville 
Distribution (WVD) analysis. This analysis was first introdueed by Wigner. The proeess 
was introduced separately by Ville. We use the WVD to produce the time-frequency 
analysis that enables the extraction of the LPI signal parameters. These time-frequency 
plots are also used to identify the modulation types. This is accomplished by sending a 
portion of the received signal through a bank of kernel functions. The length of the 
kernel function bank is determined by selecting the correct window length that provides 
the right balance of time-frequency resolution, but also stays within the computational 
limits of the system. The Fourier transform is taken of the kernel output to yield the 
frequency portion of the plot. 

This analysis contains “cross terms,” which are unintended constructive 
interference patterns that result from the presence of two signals. The locations of the 
interference patterns are desirable for signal analysis and carrier frequency 
identifications. The graphics are in three dimensions (time, frequency, and magnitude) or 
they can be displayed in a two dimensional format of time versus frequency with a color 
scheme used to show relative magnitude differences. With these plots, the modulation 
period, center frequency, and time/frequency resolutions are clearly identifiable [1]. 

2. LPI Waveform Parameter Extraction Capability 

a. Orthogonal Polyphase Sequence of 40 Sub-codes 

The WVD time-frequency image of the OPS 40 signal code set #1 with a 
center frequency of 1 kHz and bandwidth of 1 kHz is shown in Figure 58 (see Table 2). 
The WVD marginal frequency domain image Hz is shown in Figure 59. In the marginal 
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frequency domain image, the center frequency is readily apparent. In addition, the 
bandwidth of the signal can also be estimated. 


Contour of WD for O 40 4 4 C1-1-7-40-1-s 



Figure 58. Time Frequeney WD Plot OPS 40 Code 1. 
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Figure 59. Marginal Frequency Domain Plot OPS 40 Code 1. 
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b. Orthogonal Polyphase Sequence of 128 Sub-codes 


The WVD time-frequency image of the OPS 128 signal code set #1 with a 
center frequency of 1 kHz and bandwidth of 1 kHz is shown in Figure 60 (see Table 4). 
The WVD marginal frequency domain image is shown in Figure 61. In the marginal 
frequency domain image, the center frequency is readily apparent. In addition, the 
bandwidth of the signal can also be estimated. 


CMoiromD ta-o lasa 4 



Figure 60. Time Frequency WD Plot OPS 128 Code 1. 
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Mesh ofWD for 0 128 3 4 C1-1-7-128-1-S 
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Figure 61. Marginal Frequency Domain Plot OPS 128 Code 1. 

c. Optimized Orthogonal Polyphase Sequence 

The WVD time-frequency image of the Optimized OPS 40 signal code set 
#1 with a center frequency of 1 kHz and bandwidth of 1 kHz is shown in Figure 62 (see 
Table 6). There is an identifiable pattern on the slope of half of the bandwidth (BW) over 
half of the code period (T). The WVD marginal frequency domain image is shown in 
Figure 63. 
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Figure 62. Time Frequency WD Plot Optimized OPS 40 Code 1. 



Figure 63. 


Marginal Frequency Domain Plot Optimized OPS 40 Code 1. 
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d. Discrete Frequency Coding 

The WVD time-frequency image of the DFC 32 signal code set #1 with 
Np = 32, base frequency is / = IkHz, = 100 kHz, and tp = 0.001 s is shown in Figure 

64 (note the reversal of the time axis). The WVD marginal frequency domain image of 
the DFC 32 signal code set #1 is shown in Figure 65. The WVD images yield few details 
about the DFC 32 signal. The frequencies can be identified (see Table 8, Code 1) up 
until the point where the actual frequency values are masked by the cross terms. 


Cfirfni If WDfer DFC»9-1-100-1000-8 



Figure 64. Time Frequency WD Plot DFC 32 Code 1. 

Note the first five frequency multiplier terms (from Table 8) in reverse order. 
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Mesh of WD for DFC-32-1-100-1000-S 
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Figure 65. Marginal Frequency Domain Plot DFC 32 Code 1. 


3. Conclusions 

The WVD process appears to provide some signal extraction capabilities against 
the OPS signals. In all OPS simulations center frequency was readily apparent. The 
DFC signal’s parameters were nearly impossible to extract with the WVD process. Due 
to the high number of frequencies in the sequence the WVD process was able to identify 
patterns created by the waveform, but useful parameter extraction could not be 
performed. 

B. QUADRATURE MIRROR FILTER BANK ANALYSIS 
1. Introduction 

The second waveform technique used on the LPI waveforms was the Quadrature 
Mirror Filter Bank (QMFB) analysis. This analysis uses a tree of wavelet filters to 
dissect the incoming LPI signal into a series of time-frequency layers. The lower layers 

have small time resolution, but large frequency resolution. The upper layers have the 
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opposite characteristics with large time resolution, but small frequency resolution. The 
simultaneous analysis of multiple layers allows for a very accurate depiction of the 
characteristics of the LPI signal. The decomposition of the received signal is best looked 
at as a myriad of bandpass filters called wavelets that are used to break up the signal into 
smaller pieces that are better for analysis. This is accomplished by using a standard 
wavelet as the basis function and then translating, dilating, or contracting the wavelet as 
necessary to achieve the desired time and frequency resolution. The signal parameters 
are graphically displayed in time-frequency plots that have color schemes to identify 
centers of concentrated signal energy. The QMFB has the advantage over WVD and 
Cyclostationary Spectral Analysis (CSA) because it has no cross terms that adds noise to 
the signal and could possibly obscure extraction of its parameters. The output graphics 
are used to identify modulation period, center frequency, signal bandwidth, code period, 
and in some cases the phase shift [1]. 

2. LPI Waveform Parameter Extraction Capability 

a. Orthogonal Polyphase Sequence of 40 Sub-codes 

The fourth layer of the QFMB time-frequency plane of the OPS 40 signal 
code set #1 with a center frequency of 1 kHz and bandwidth of 1 kHz is shown in Figure 
66 (see Table 2). Five periods are captured. The repetitive nature of the polyphase 
sequence and the center frequency are obvious as shown in Figure 66. A close up of the 
fourth layer of the QFMB time-frequency plane of the same OPS 40 signal is shown in 
Figure 67. The code period is readily apparent in this view. 
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Figure 66. QMFB Layer 4 of OPS 40 Code 1 showing all five eode periods. 



1 0.9 
0.8 
0.7 
0.6 
0.5 
'0.4 

1 0.3 
0.2 

0.1 



Figure 67. 


QMFB Layer 2 of OPS 40 Code 1 with eode period T = 40 ms. 
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b. Orthogonal Polyphase Sequence of 128 Sub-codes 


The fourth layer of the QFMB time-frequency plane of the OPS 128 signal 
code set #1 with a center frequency of 1 kHz and bandwidth of 1 kHz is shown in Figure 
68 (see Table 4). Five periods are captured. The repetitive nature of the polyphase 
sequence and the center frequency are obvious in Figure 68. A close up of the fourth 
layer of the QFMB time-frequency plane of the same OPS 128 signal is shown in Figure 
69. The code period is readily apparent in this view. 
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Figure 68. QMFB Layer 4 of OPS 128 Code 1 showing all five eode periods. 
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Figure 69. 


QMFB Layer 4 of OPS 128 Code 1 showing code period T = 128 ms. 


c. Optimized Orthogonal Polyphase Sequence 

The fourth layer of the QFMB time-frequency plane of the Optimized 
OPS 40 signal code set #1 with a center frequency of 1 kHz and bandwidth of 1 kHz is 
shown in Figure 70 (see Table 6). Five periods are captured. The repetitive nature of the 
polyphase sequence and the center frequency are obvious in Figure 70. A close up of the 
fourth layer of the QFMB time-frequency plane of the same Optimized OPS 40 signal is 
shown in Figure 71. The code period is readily apparent in this view. 
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Figure 71. 
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QMFB Layer 4 of Optimized OPS 40 Code 1 all five code periods. 



QMFB Layer 2 of Optimized OPS 40 Code 1 code period T = 40 ms. 
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d. Discrete Frequency Coding 

The fifth layer of the QFMB time-frequency plane of the DFC 32 signal 
code set #1 is shown in Figure 72 (see Table 8). A close up of the fifth layer of the 
QFMB time-frequency plane of the same DFC 32 signal is shown in Figure 73. The 
frequency duration is discemable in this view and the frequency values proportional to 
the original code set can be determined. QMFB presents a clearer view of frequency 
duration and frequency value than WVD. 
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Figure 72. QMFB Layer 5 of DFC 32 Code 1. 
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Figure 73. The First Period of QMFB Layer 5 of DFC 32 Code 1. 


3. Conclusions 

QMFB showed some potential for parameter extraction of the OPS code. The 
code period was identifiable utilizing QMFB under no noise simulations as well as an 
estimate of center frequency. The QMFB process performed on the DFC signal revealed 
that a frequency hopper was being utilized, and the extraction of the sequence is made 
possible by this method. In addition, with a close in view the frequency duration could 
be determined. 
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C. CYCLOSTATIONARY SPECTRAL ANALYSIS 

1. Introduction 

The third waveform analysis teehnique used on the LPI signals was the 
Cyclostationary Spectral Analysis (CSA). The main premise of CSA is the analysis of 
the received signal using a different less traditional approach. Both WVD and QMFB 
viewed parameter extraction from a time-frequency point of view. CSA uses a frequency 
versus cycle frequency format in order to provide an additional capability for LPI signal 
analysis. This bi-frequency analysis is obtained through using Fast Fourier Transforms 
(FFT) to manipulate the ACF into a cyclic autocorrelation function (CACF). The Fourier 
transform is applied to the CACF to produce the CSA. The CSA can be found with either 
a time-smoothing FFT Accumulation Method (FAM) or a Direct Frequency-Smoothing 
Method (DFSM). The FAM uses many smaller calculations that produce a more timely 
analysis of the signal. DFSM is more computationally intensive and relies on spectral- 
correlation operations. The plots generated are used to identify center frequency, signal 
bandwidth, and code rate [1]. For consistency all CSA was conducted selecting a 
frequency resolution of 128 Hz and a df/da (M) of 2. 

2. LPI Waveform Parameter Extraction Capability 

a. Orthogonal Polyphase Sequence of 40 subcodes 

The frequency smoothing spectral correlation density (SCD) of the OPS 
40 signal code set #1 with a center frequency of 1 kHz and bandwidth of 1 kHz is shown 
in Figure 74. A closer look, shown in Figure 75, at a single modulation quadrant reveals 
additional signal structure. The time smoothing spectral correlation density (SCD) of the 
OPS 40 signal code set #1 with a center frequency of 1 kHz and bandwidth of 1 kHz is 
shown in Figure 76. A closer look, shown in Figure 77, at a single modulation quadrant 
reveals additional signal structure. 
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Frequency Smoothing SCD O 40 4 4 C1-1-7-40-1-s, df = 128, N = 128 



Figure 74. DFSM of OPS 40 Code 1 with T = 40 ms and B = 1000 Hz. 


Frequency Smoottirv SCD 0 40 4 4 C1-1-7-40-1-S. df = 128, N = 128 



Figure 75. DFSM of OPS 40 Code 1 with T = 40 ms and B = 1000 Hz Individual 

Modulation Quadrant. 
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Time Smoothing SCD O 40 4 4 C1-1-7-40-1-s, df = 128, N = 128 



Figure 76. FAM of OPS 40 Code 1 with T = 40 ms and B = 1000 Hz. 


Time Smooltirxi SCD 0 40 4 4 C1-1-7-40-1-6. df = 128, N = 128 



Figure 77. FAM of OPS 40 Code 1 with T = 40 ms and B = 1000 Hz Individual 

Modulation Quadrant. 
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b. Orthogonal Polyphase Sequence of 128 Sub-codes 


The frequency smoothing spectral correlation density (SCD) of the OPS 
128 signal code set #1 with a center frequency of 1 kHz and bandwidth of 1 kHz is shown 
in Figure 78. The center frequency can be observed in Figure 78. A closer look, shown 
in Figure 79, at a single modulation quadrant reveals additional signal structure. The 
time smoothing spectral correlation density (SCD) of the OPS 128 signal code set #1 with 
a center frequency of 1 kHz and bandwidth of 1 kHz is shown in Figure 80. A closer 
look, shown in Figure 81, at a single modulation quadrant reveals additional signal 
structure. 


Frequency Smoothing SCD O 128 3 4 Cl-1-7-128-1-s, df = 128, N = 128 
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Figure 78. DFSM of OPS 128 Code 1 with T = 128 ms and B = 1000 Hz. 
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Frequency SrnootMr« SCO O 128 34 C1-1-7-128-1-S. df = 128, N = 128 



Figure 79. DFSM of OPS 128 Code 1 with T = 128 ms and B = 1000 Hz Individual 

Modulation Quadrant. 


Time Smoothing SCD O 128 3 4 C1-1-7-128-1-s, df = 128, N = 128 



Figure 80. FAM of OPS 128 Code 1 with T = 128 ms and B = 1000 Hz. 
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This SmootMng SCO O 128 3 4 C1-1-7-128-1-8. df = 128, N = 128 



Figure 81. FAM of OPS 128 Code 1 with T = 128 ms and B = 1000 Hz Individual 

Modulation Quadrant. 


c. Optimized Orthogonal Polyphase Sequence 

The frequency smoothing spectral correlation density (SCD) of the 
Optimized OPS 40 signal code set #1 with a center frequency of 1 kHz and bandwidth of 
1 kHz is shown in Figure 82. The center frequency can be observed in Figure 82. A 
closer look, shown in Figure 83, at a single modulation quadrant reveals additional signal 
structure. The time smoothing spectral correlation density (SCD) of the Optimized OPS 
40 signal code set #1 with a center frequency of 1 kHz and bandwidth of 1 kHz is shown 
in Figure 84. A closer look at a single modulation quadrant as shown in Figure 85 
reveals additional signal structure. The energy centers present indicate a distinct linear 
modulation that is present in orthogonal waveforms. 
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Figure 82. DFSM of Optimized OPS 40 Code 1 with T = 40 ms and B = 1000 Hz. 
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Figure 83. DFSM of Optimized OPS 40 Code 1 with T = 40 ms and B = 1000 Hz 

Individual Modulation Quadrant. 
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Figure 84. 


Figure 85. 


Time Smoothing SCD KO 40 3 4 CM-7-40-1-s, df = 128, N = 128 



FAM of Optimized OPS 40 Code 1 with T = 40 ms and B = 1000 Hz. 


TimeSmoolMngSCO KO4034C1-1-7-40-1-S,df = 12B, N= 128 



FAM of Optimized OPS 40 Code 1 with T = 40 ms and B = 1000 Hz 
Individual Modulation Quadrant. 
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d. 


Discrete Frequency Coding 


The frequency smoothing spectral correlation density (SCD) of the DFC 
32 signal code set #1 is shown in Figure 86. A closer look, shown in Figure 87, at a 
single modulation quadrant reveals additional signal structure. The time smoothing 
spectral correlation density (SCD) of the DFC 32 signal code set #1 is shown in Figure 
88. A closer look, shown in Figure 89, at a single modulation quadrant reveals additional 
signal structure. For CSA on the DFC32 signal the base frequency was reduced to 50 Hz 
and the sampling frequency was reduced to 5 kHz. 


Frequency Smoothing SCD DFC-32-1-5-1000-S, df = 128, N = 128 



Figure 86. DFSM of DFC 32 Code 1 with / = 1000 Hz. 
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Frequency Smoothing SCD DFC-32-1-5-1000-S, df = 128, N = 128 
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Figure 87. DFSM of DFC 32 Code 1 with / = 1000 Hz Individual Modulation 

Quadrant. 


Time Smoothing SCD DFC-32-1-5-1000-S, df = 128, N = 128 



Figure 88. FAM of DFC 32 Code 1 with / = 1000 Hz. 
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Time Smoothing SCD DFC-32-1-5-1000-S, df = 128, N = 128 



Figure 89. FAM of DFC 32 Code 1 with / = 1000 Hz Individual 

Modulation Quadrant. 


3. Conclusions 

Cyclostationary spectral analysis provides the most insight to center frequency 
and bandwidth measurements of the LPI signals. Each pattern was anchored by a 
multiple of the center frequency and an individual section of the pattern was twice the 
width of the bandwidth. Within each pattern, the code rate and code period can also be 
measured to obtain a better idea of the coding being used and the origin of the waveform. 
The polyphase signals with > 40 complicate the extraction of detail down to the sub¬ 
code level. 
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V. UAV SWARM THEORY 


A. UAV SWARMS 

The concept of a swarm of individuals acting in conjunction with one another to 
accomplish a greater task is not a new idea. Entire colonies of animals, insects, and 
organisms depend on the foundations of a swarm to survive. Take ants fighting off a 
predator for example. A predator many times larger and stronger than an individual ant is 
easily subdued by the coordinated efforts of a swarm of ants too numerous for the single 
predator to fend off. A pride of lions hunts in a pack or “swarm.” They attack animals 
much more powerful and sometimes faster and more capable, but with their combined 
efforts, they capture their prey. From a different viewpoint, the lions’ prey by themselves 
does not have a chance to survive against the pack of lions. However, as a herd or 
“swarm”, the herd always survives because of their numbers and coordinated escape. 
From nature it can be observed that three main pillars support the swarm concept: 
survivability of the swarm, the coordination of the swarm through internal 
communication, and enhanced capabilities of the swarm versus individual capabilities. 

The UAV has significant advantages over a traditional surveillance aircraft. First, 
their name implies that there is no risk of loss of human life when a UAV is sent to a 
battlefield to perform its mission. The price tag on a UAV is also significantly less than a 
fighter aircraft, which allows the UAV to enter situations that pose more danger because 
the risk assessment is lessened. Second, a UAV can operate in an area that is not suitable 
for human life due to nuclear, chemical, or biological weapons that have been detonated 
there. This provides for battle damage assessment and other missions of that kind. 
Finally, the radar cross section (RCS) of a UAV is quite small compared to a fighter 
aircraft. This means the UAV can infiltrate enemy defenses much further than a fighter 
and perform missions closer to the enemy [20]. 
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B. SWARMING MECHANISM 


The swarming mechanism is based on the ant-inspired model of foraging for food, 
which is also displayed in many other biological creatures. The swarm is comprised of 
individual entities or small units that are in search of a target of interest. The search for 
such a target is coordinated with the use of chemical substances or pheromones. 
Pheromone trails are used to guide other ants to the spot where a target of interest has 
been spotted. Pheromones are continuously laid as the ant searches for a target. As the 
ant moves further from the colony, the strength of this scent of its trail diminishes. Once 
the target (use food for instance) has been spotted, the ant takes some of the food and 
brings it back to the colony. The return trip reinforces the pheromone trail thus 
confirming the validity of the target. As more and more ants follow the trail to the food, 
the pheromone trails are strengthened and it is easier for successive ants to gather food 
for the colony [11]. 

Termites are another creature that can build elaborate mud arches using the 
principle of pheromones. The termite senses the pheromone, which tells him to gather up 
dirt in its mouth and deposit it where the scent of the pheromone is the strongest. The 
termite gathers dirt and moves to the location of the strongest scent and deposits the mud 
there. As this process is repeated by a number of termites all leaving their own 
contribution of pheromone at the scene, the autocatalytic effect takes over. This positive 
feedback cycle reinforces the strength of the pheromone scent and lures more and more 
termites to the scene to contribute to the building of the arch, which in turn brings the 
arch closer to completion. Once the arch has been completed, there will be no more 
termites following the trail to build the arch, so the pheromones will wear out and the 
scent will vanish [13]. 

Fixed action patterns and instinctive behaviors can be explained using the concept 
of pheromones. These chemicals leave behind odors that, in turn, signal or alert other 
ants. They are drawn to the possible location based on what signal that particular 
pheromone sends. Various pheromones are used to trigger different simplistic responses 
from ants. Pheromones are used to adapt the behavior of the individuals into group 
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behavior that is more useful and stronger than that of the individual ant. Hence, 
pheromones are a means of communication between ants that allows them to coordinate 
their efforts to enable their group to accomplish a goal. Individually, the ants are not a 
force to be reckoned with in the least, but coordinated together to form a swarm or group; 
the ants can accomplish extreme feats [11]. 

The main idea in swarm intelligence is the accomplishment of a grand feat 
without centralized control, but by using simple tasks that sum together to realize the 
whole goal. This decentralized control is brought about by stigmergy, which is 
communication by altering the state of the environment in such a way that will affect the 
behaviors of others for whom the environment is a stimulus. There are two different 
categories of stigmergy, cue-based and sign-based. Cue-based stigmergy defines the way 
a certain event or action is a tell tale indication that the behavior of the individual should 
change in order to conform to the idea of the group. Sign-based stigmergy refers to the 
signal that is sent to the other individuals, possibly by way of pheromone trail, to change 
to their actions. The arch building termites’ example includes both cue-based and sign- 
based stigmergy. The cue-based stigmergy is seen as the tower building process. A 
termite will see the tower building process beginning and know to gather dirt in its mouth 
to contribute to the building efforts. The sign-based stigmergy is the signal of the 
pheromone of where to place the dirt from the mouth of the termite. Both stigmergetic 
effects are geared towards the same goal and continue until the task has been completed 
[13]. 


C. UAV SWARM CONCEPT 

UAV have become an integral part of the modern battlefield. They have a variety 
of missions ranging from Intelligence, Surveillance, and Reconnaissance (ISR) missions, 
to directed fire support of ground troops, to unmanned strike missions. The versatility of 
the UAV, its relatively low cost as compared to a fighter plane, and its ability to remain 
hidden from enemy radar (low radar cross-section) make it a very useful platform. 
However, by itself, the UAV can be overcome by enemy forces. The UAV has little to 
no armor or defensive components that will ensure it survivability in the case of attack. 
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To combat this, the idea of the UAV swarm eame about to ensure the survivability of the 
UAV in order to accomplish the mission. The inereases produetivity and output of the 
UAV swarm also aided in the development of UAV swarm theory [19]. 

UAV swarm theory is explained as offensive action generated in pulses by highly 
dispersed forees that do not employ traditional hierarchical command control structures. 
The 9/11 attaeks are an example of this swarm doctrine. There are many different 
separate entities that are eapable of performing the mission by themselves. The 
individual UAV are given simple instructions such as, “perform reeonnaissance on an 
enemy military eompound.” Each UAV can accomplish this mission on its own; 
however, the more UAV engaged in the mission, the faster it is eompleted and the better 
the finished product. 

Historical analyses have shown that the Objective Force must possess superiority 
in the following three battlefield functional capabilities in order to defeat future threats: 
situational awareness, stand-off lethality, and survivability by agility. Situational 
awareness is derived from the timely processing and fusion of remote sensor data. Stand¬ 
off lethality is highly dependent on an effective merge of sensor grids and shooter grids 
as targeting data will be eontained on both grids and is neeessary for an aeeurate and 
timely shot. The survivability of a team of UAV is neeessary for the post-modem 
battlefield. Knowledge of when and where to move will minimize vulnerability to enemy 
defense networks and maximize the capabilities of the autonomous platform [12]. 

The core idea of the UAV swarm is eentered on the ability of the swarm to be 
able to communieate within itself in order to eoordinate their efforts to maximize 
efficieney and minimize their vulnerability to attack. The possible solution to this is the 
formation of a “gossip network” between UAV that is used to transfer information from 
one UAV to another. This wireless network will be open to any UAV joining the mission 
and will be the medium through whieh the Digital Information Pheromones (DIP) will be 
transferred to alert the other UAV in the swarm about the presenee of a target or a threat. 
Mueh like the ants and termites mentioned above, the pheromones will be digitized and 
sent through the gossip network and serve the same purpose of alerting other UAV to a 


discovery made by a UAV in the swarm [19]. 
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The UAV that is uploading important DIP on a possible target location to the 
pheromone landscape broadcasts a message that alerts other UAV. This message is 
received by the other UAV in range and directs them to subscribe to the gossip network 
alerting them to the presence of a new post. The other UAV in range of the initiator 
UAV will respond with a message acknowledging they have received the message and 
will establish the wireless communications channel. Each message will contain the 
location and message time of the initiator UAV. The probability of passing on gossip 
messages and posts decreases exponentially as distance from the initiator UAV increases. 
This can be counteracted by the UAV that receive the initiator message relaying that 
message to other UAV that might be further away from the initiator UAV. This network 
architecture guarantees the survivability of the gossip network and allows for any of the 
UAV in the swarm to be a hub, relay, or member in order to complete the mission [19]. 

D. ADVANTAGES OF THE UAV SWARM 

The advantages of the UAV swarm are precipitated in the interaction within the 
swarm and the speed with which the swarm of UAV can coordinate their efforts to 
accomplish their mission. The use of Internet Protocol to convey information over a 
network has inherent advantages that the UAV gossip network will build on [12]. The 
use of flags or indicators as to the purpose or content of a packet is used often in Internet 
Protocol to alert different machines as to the nature of the information they are receiving. 
The flags that would be most useful to the UAV swarm would be: Deploy (in transit to 
search region). Exploring (searching region for potential targets). Following (leaves 
search and traverses path indicated by digital pheromone path to potential target). 
Monitoring Target (flying close enough to the target to survey it), and Return to Base 
(UAV has remaining fuel lower level interlock that forces it to disregard further gossip 
network flags and proceed back to base). Much like the biological creature pheromones, 
these five flags would be necessary to direct the simple actions of the UAV in order to 
contribute to the swarm that is trying to accomplish it mission 13]. 

Two of the flags that can be further manipulated to yield more accurate direction 
to the UAV are the Following and Monitor Target flags. The Following flag is used at 
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the onset of target identifieation by one of the UAV in the swarm. The UAV issues the 
flag in a message, but there are eonditions on the flag released. These eonditions refer to 
the relative eertainty with whieh the initiator UAV is of finding a potential target. If the 
UAV has deteeted an enemy radar signal of interest in the threat band and held the signal 
for enough time to analyze and extract its parameters, then the Follow flag will have a 
high importance attached to it. This signifies to the other UAV that they should cease 
their Exploring flag and follow the DIP to the location of the initiator UAV to begin the 
Monitor Target flag. The other UAV send their messages back to the initiator UAV then 
head to the first waypoint indicated by the DIP. From here, the DIP will lead the UAV to 
the location of the first possible target contact and begin their Monitor Target sequence. 
Conversely, if there was only a brief hit of a radar signal that was unable to be analyzed 
and therefore only a sporadic indication of a possible target, there will be a low 
importance flag placed on an uncertain target sighting. Only the closest UAV might 
leave its search and follow the DIP to the target sighting [13]. 

The Monitor Target flag is a short duration flag that is used only when a target has 
been found and data collection is required to fix its position and record its location for 
future analysis. The Monitor Target flag will be switched on once the UAV arrives at the 
end of the DIP trail and acquires the possible target’s signal. The UAV will then take up 
a monitoring path that resembles a 90-degree arc around the potential location of the 
target. The UAV will travel along this arc ensuring it keeps constant contact with the 
target signal, but remaining far enough away from the target as to keep from arousing 
suspicion. The UAV will collect data and attempt to extract the parameters of the target 
signal using its onboard software. Once the UAV has completed the arc, it will have an 
initial bearing to the signal and a cross bearing from the final bearing to the target. The 
place where these bearings cross will be recorded by the UAV and will be stored as the 
location of the target emitting the signal. This data will be uploaded to the gossip 
network and shared with other UAV to find the most probable location of the target [12]. 
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E. CHALLENGES TO THE SWARM CONCEPT 


Despite the gains in surveillance and reco nn aissance that a swarm of UAV has 
over a single platform performing a similar mission, there are some challenges that must 
be addressed. First is the lack of a physical medium on which to leave a pheromone trail. 
The first solution to this challenge is somewhat complicated, but given enough 
development of the idea could provide a nearly untraceable means of communication 
between UAV. When a UAV wishes to make its message known to the rest of the group, 
it releases a chemical into the air. This chemical will have much of the same properties 
as air (density, molecular composition, and temperature) to ensure it does not float away 
or sink to the ground. Each UAV will have onboard a miniature gas chromatographer 
that can separate the atomic elements in the air and analyze them. When another UAV 
detects the chemical that signifies something to the swarm of UAV, such as target has 
been spotted, it will follow the trail of the chemical until it rendezvous the initiator UAV 
at the target sighting location. This method of communication is quite primitive and 
requires further development before bringing it into service, but it most closely resembles 
that of the biological pheromone trail [13]. 

The other possible solution to the physical pheromone medium challenge has 
been introduced above and that is the idea of a gossip network. This network would be 
localized to the UAV in the swarm and would be used to transfer information and alerts 
to the other UAV in the swarm. Using Internet protocol and its use of “flags”, the UAV 
will communicate and share information on the network. Each UAV will be programmed 
to periodically subscribe to the gossip network to receive updated tasking and changes to 
the mission. In addition, target sightings, locations, and signal parameters of interest will 
be available for UAV to peruse when they get a chance. Using this network, the UAV 
will also alert one another to the possible sighting of a target and using DIP they can 
guide each other on to the trail of the target’s signal. This coordination is what gives the 
swarm of UAV the advantage over one UAV by itself [12]. 

A second challenge that will have to be addressed is the possibility that a number 
of the targets being observed by the UAV swarm may not be fixed. Normally, ant targets 
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are fixed, but UAV targets ean be mobile. Therefore, it has been suggested that an anti¬ 
pheromone be used to update the loeation of mobile and dynamie targets. This means 
that in addition to the DIP, there will be a eancellation DIP that will negate the previous 
pheromone. On the other hand, if the trail to a target changes, ants continue the trail of 
DIP to the new location of the target in order to regain access for the rest of the swarm. 
This continued trail if often inefficient and many times more lengthy that one continuous 
trail. This aspect of the discussion is critical due to the limited fuel capacities of UAV, 
the concern for weight restrictions, and the reluctance to take up cargo space with extra 
fuel. Anti-DIP is a much more direct and efficient way to make a path to a target, but 
continuing the DIP trail from where it left off to the new location of the target is faster 
and can be updated closer to real time [12]. 

F. CONCLUSIONS 

The idea of taking a group of advanced unmanned machines and coordinating 
them to work together like a Navy SEAL Team or a Marine Scout Sniper is a very 
serious undertaking. With a human unit, there are many concerns such as survivability, 
sustainability, security, shelf-life of the asset, and most important of all the loss of human 
life. Many of these concerns are alleviated by using a swarm UAV to accomplish the 
same mission, but with real time analysis of the intercepted signal, a more efficient and 
timely acquisition of the target, and above all the absence of the risk of human life. The 
technology for DIP and gossip networks are progressing and will aid the military in their 
future ISR missions. 
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VI. UAV SWARM JAVA PROGRAM 


A. PROGRAM BACKGROUND 

1. Program Development 

The Java program that came before this program was originally created by Chin 
Lua of the Computer Science and Operations Research Department at North Dakota State 
University with the intent on employing a swarm of UAV to partake in a synchronized 
attack against a target of opportunity [21], The UAV swarm is sent out with the mission 
to destroy a certain target (enemy radar or other transmitter). Once one of the UAV 
detects the target, all the other UAV are notified and a surveillance mission is conducted. 
Once, the target has been adequately confirmed to be the target of interest, the UAV take 
up attack positions and wait for the countdown to strike. Each UAV carries a small 
payload and all the munitions are dropped simultaneously on the target at the 
commencement of the attack. 

This Java program has a slightly different mission and therefore has different 
goals and an entirely different end state. The swarm of UAV is not used for combat or 
strike missions, but as an ISR mission capable group. The swarm of UAV will be 
uniquely outfitted with the signal analysis tools discussed above in order to perform ISR 
missions on all LPI systems and radar signals. The threat of radar systems and netted 
radar using LPI waveforms a great one for strike aircraft whose radar cross-section is 
enormous compared to that of a UAV. The Java program shows a simulated ISR mission 
of a UAV swarm that has been deployed to hunt down an element of an Integrated Air 
Defense System (IADS). These radar systems could be monostatic or bistatic depending 
on the technology and the level of integration of the radar network. 

The UAV will set out from their base rendezvous point and search an area for LPI 
signals that correspond to known threat bands of enemy LPI radar. Once found, the UAV 
detecting the enemy signal will broadcast the detection to the other UAV in the area to 
assist in data collection and target confirmation. The UAV will take up stations called 
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orbit stations in an arc at the maximum detection range from the target for survivability 
concerns. Once the UAV have subscribed to the gossip network and posted their plotted 
position of the enemy radar target, they will move on. This will continue until: The 
mission is complete and the UAV are instructed to return home; the low-fuel flag is 
detected by the swarm and they are ordered to return home; or they are engaged by 
enemy fire and the mission is aborted. 

2. UAV Response Actions 

The UAV swarm will employ three basic actions during its mission. These 
actions are implemented in the hardware of each UAV, so communication within the 
swarm for contact avoidance is not necessary. The first of these is the “avoid” maneuver. 
It entails a UAV detecting two other UAV or any obstacles in its path. The UAV will 
perform a simple vector analysis of the situation and act accordingly to avoid the two 
obstacles. The UAV will use the bearings to each of the two obstacles summed together 
to form the resulting course for contact avoidance. The UAV will continue to move 
along the new vector and turn clockwise or counterclockwise to avoid any obstacles. 
This iteration of vector analysis is computationally simplistic in nature and once the UAV 
is clear of the obstacles, it will return to its original mission of searching for a target or 
observing a target. 

The second simple UAV action is the “attract” maneuver. This movement is in 
response to a UAV receiving confirmation of a target sighting. The UAV will advance 
towards the target until it also detects the target. The UAV plots a set of waypoints in a 
circular fashion around the target. Next, the UAV selects the waypoint most directly 
between the target and itself and promotes this waypoint to make it the station point. 
Finally, the UAV moves to the station point in the set to allow it to properly position 
itself for data collection and to complete its ISR mission. 

The third action the UAV performs is the “orbif’ maneuver. This movement 
consists of the UAV finding its station point around which to set up a circular track and 
then traversing that track or orbit while it completes its mission. The UAV plots a set of 
waypoints in a circular path around the station point in one quadrant about the target. 
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The UAV then follows these plotted waypoints while it gathers intelligenee on the target. 
This intelligence will be used to properly identify the target once the mission is 
completed and most importantly to fix its position for mapping and possible targeting for 
destruction at a later time. The UAV continues to move in a counterclockwise pattern 
about the reference point at the maximum detection range of the UAV to the target. Once 
the requisite amount of intelligence has been gathered, the UAV moves on to search for 
another target according to its mission parameters [15]. 


In Figure 90 shown below, the movements and actions of the individual UAV are 
conceptualized to provide a basis for their actions in the Java Program. The UAV control 
structure employs a decentralized architecture that is the basis for their autonomous 
mission capabilities. The UAV respond to their external sensors and inputs from the 
gossip network. These factors determine the next course of action for the UAV. As part 
of a UAV swarm, much like the insects addressed above, the initiation of action by one 
UAV is propagated to the other UAV in order to accomplish the mission through 
coordinated action. 



Avoid 


Artack 

(Atti'actioii) 

Orbit Station 
(Orbit) 

Orbit Target 
(Orbit) 

Search 

(Atti^actioii) 



Figure 90. UAV Control Structure in Java Program. 
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3. 


The Gossip Network Function In The UAV Swarm 


The gossip network is essential for the coordination of the UAV in accomplishing 
their mission. In Figure 91 below, the gossip network traffic is depicted. The UAV 1 is 
the UAV that in both cases detects the target first. The target detection is denoted by the 
step increase in awareness of the UAV. This awareness includes the transition from 
“search” maneuvers to “attract” maneuvers. Once the UAV has calculated the station 
points and determined which station point to observe the target from, it transitions to 
“orbit” maneuvers. The sustained value of 0.5 higher than the base value of one for UAV 
1 shows a heightened sense of awareness during which time the UAV is observing the 
target, recording bearing and range information, and passing this data to the other UAV 
over the gossip network. Once the UAV has finished the requisite actions for target 
observance, it disengages from the target as shown by the step decrease in awareness to a 
value of 0.5 lower than the base value. During this time, the UAV is turning to a course 
that takes it safely away from the target and reassuming its “search” maneuvers. 


Gossip Network Traffic 



Figure 91. Gossip Network traffic from three UAV in their search of two targets. 
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The other two UAV show the results of the message traffie taking place on the 
gossip network. UAV 2 is the first to receive the target detected broadcast message from 
UAV 1. The ramp up in awareness for UAV 2 over the course of five time slots shows 
its transition both in position relative to the target and behavior state. UAV 2 is not so 
close that it detects the target when UAV 1 does, so it must travel to the waypoint 
broadcast by UAV 1 once target detection occurs. When UAV 2 reaches a value of 0.5 
higher than its base value, it has reached the target detection waypoint and is now 
performing “attract” maneuvers in order to get into position to correctly observe the 
target. Once it has reached its station point, it begins to orbit that point and take data on 
the target by performing “orbit” maneuvers. Upon completion of data collection, UAV 
follows suit with UAV 1 and returns to its “search” maneuvers to look for the second 
target. 

UAV 3 is the furthest away from the target upon initial detection. UAV 3 does 
not hear the target detected message broadcast from UAV 1, so that is why its ramp up in 
awareness occurs one time slot after UAV 2 begins to ramp up awareness. This is 
because as part of the gossip network, UAV 2 rebroadcasts the target detected message as 
soon as it is received. UAV 3 is in range of UAV 2, so it hears the message and is thus 
delayed in reaching the target. In a similar fashion as the other UAV, UAV 3 transitions 
from “search” maneuvers to “attract” maneuvers upon reaching the target detection 
waypoint. Once its station point is calculated, it travels to the station point and begins to 
orbit the station point and perform orbit maneuvers. After sufficient data has been 
collected and transmitted through the gossip network, UAV 3 disengages from the target 
and falls back to its “search” maneuvers in order to find the next target. 

B. THE UAV SWARM JAVA PROGRAM 

1. The Java Program Equivalence 

The actions and maneuvers of the UAV swarm are based in part on the 
capabilities of the Compact LCR-100 Electronic Intelligence (ELINT) receiver that can 
be carried onboard a UAV with little restriction on overall payload and endurance. This 
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is an example of the type of sensor needed to enable the performance of the UAV swarm 
in a similar manner to the Java program. This receiver weighs about a pound and 
occupies only 12 cubic inches of payload space. Its sensitivity is approximately -105 dB 
and covers a frequency range of 0.5 - 18 GHz. These statistics are a minimum 
requirement to ensure the effectiveness and survivability of the UAV Swarm. 

All the measurements and distances in the UAV Swarm Java program are 
referenced to pixels on a computer screen. The target mark that the UAV are trying to 
detect is a diameter of 10 pixels. To make a more useful parallel, it is assumed that a 
radar site would be 150 feet in diameter; hence, the ratio of computer pixels to feet is 
1:15. The vision range of a UAV searching for a target is 400 pixels or 6,000 feet. This 
is just over a mile, which is typical when searching for a LPI emitter. This means that the 
swarm of UAV must comb their search terrain and come to within a mile of the LPI 
emitter in order for the ELINT receiver to detect the LPI signal waveform and extract its 
parameters using the methods discussed above. This range is very close to the LPI 
emitter, but with such a small RCS, the UAV has a low probability of detection. 

The other parameters of the UAV Swarm such as maximum speed and maximum 
turn are arbitrary constraints for a generic UAV. Some newer UAV can go faster and 
have a smaller turning radius, but these limits were chosen to maximize on station time 
and minimize fuel consumption. The repulsion range variable refers to the minimum 
distance required between one UAV and another to prevent collisions and provide for 
adequate maneuver room to avoid other contacts. A UAV model with more sophisticated 
avoidance algorithms or faster processors could reduce the repulsion range, but this value 
was chosen for the generic UAV model. 

The on station time required for adequate surveillance and target location 
determination is based in part on the time required to pull the signal out of the 
environment and in part on the direction finding (DF) and ranging capabilities of the 
UAV antenna supported by the ELINT receiver. The time required to fix the target’s 
position is inversely proportional to the number of UAV attempting to find the target. 
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This means that as the number of UAV increases, the time required to fix the target’s 
position goes down. This is another reason that a swarm of UAV is very effective in both 
finding and pinpointing a target’s position. 

The time required to fix a target’s position is largely dependent on two factors: 
the bearing and range errors of the ELINT sensor and the number and location of data 
points uploaded to the gossip network, which is dependent on the number of UAV in the 
swarm. The bearing and range error of an ELINT sensor’s DF capabilities are rarely 
adjustable. The more advanced the sensor, the smaller the error. 

The UAV swarm presents a solution to this problem by providing more data 
points in a timelier manner. One UAV trying to DF a target would produce a circular 
error probable (CEP) many times larger than a swarm of five UAV. This is due to 
inability of one UAV to get data points quick enough and accurately enough. In Figure 
91 below, the relative differences between the CEP of one UAV and five UAV are 
shown. With five UAV, the CEP is constricted due to the increased number of data 
points and the locations of the UAV. Figure 92 below shows the reduction in the CEP 
from just two UAV at crossed bearings. The gray ellipses depict the CEP produced when 
one UAV attempts to DF the target. The blackened area where the two ellipses cross 
shows the reduced area of commonality between the two bearings and ranges. The black 
circle around the target (red circle inside the blackened area) shows the resulting CEP 
when the two data points are uploaded to the gossip network. The gossip network is 
instrumental in the comparing and contrasting of the data points from all the UAV in the 
swarm. The individual UAV contribute to the final target location CEP and that location 
is saved on board every UAV to ensure survivability of the mission and use in future 
missions [22]. 
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Figure 92. Illustration of UAV Swarm reduction target location CEP [22], 

2. UAV Swarm Detection Process 

The Java Program illustrates the interaction of the UAV as they traverse their 
search area in order to accomplish their mission. Each time new data is obtained by any 
of the UAV in the swarm, the information is relayed by the gossip network to each of the 
other UAV. This promotes the coordination of UAV actions so that the effective efforts 
of the UAV swarm are increased. The individual UAV can find all the targets by itself, 
but would have to remain on station for longer in order to obtain a more accurate target 
location. Figure 93 illustrates the coordination of the UAV through the gossip network in 
finding and fixing the positions of the two targets below. 
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In the simulation in Figure 93 below, three UAV are searching for two targets 
situated about ten miles apart. UAV 1 is the black track, UAV 2 is the blue track, and 
UAV 3 is the red track. The UAV start their mission at the far left vertically stacked 
using search mode to find the first target. As discussed, this is a simple mode that tells 
the UAV to continue along a predetermined path (in this case it is a saw tooth course to 
the southeast) until it detects the target. Once target detection is accomplished by UAV 
1, its detection is relayed over the gossip network by a message containing waypoint data 
and initial radar signal parameters (these are rough estimates of carrier frequency and 
bandwidth) and broadcasts that to the other UAV in range. The other UAV broadcast the 
same message including the data time group (DTG) so each UAV can determine its 
location relative to the target detection sighting. Once a UAV receives a target detection 
message from another UAV, it alters its course to align with the target detection waypoint 
and proceeds to that point. 

Once a UAV reaches the initial target detection waypoint, it proceeds to follow its 
programmed path to determine which station point it should take up in order to perform 
surveillance and data collection on the target. Each successive UAV fills in the next 
station circle over from the previous one and begins its orbit. The maneuvers required for 
a UAV to travel from its search point to its station circle are quite simple and are mostly 
comprised of avoidance maneuvers to prevent collisions between UAV. Once each UAV 
has obtained a station circle it begin to orbit and transitions to the data collecting phase of 
the target location process. 

At the onset of a UAV reaching a station and orbiting the station point, it begins 
taking data on the target signal’s waveform. The goals are two fold: attempt to extract 
the LPI signal parameters during surveillance for future use against the target weapons 
system and obtain a requisite number of DF values to include bearing and range at a 
minimum in order to provide enough data to upload an accurate value to the gossip 
network for comparison with the other UAV in the swarm. The accomplishment of the 
first task can be difficult and possibly not occur at all. Since the UAV is gathering the 
data, it can be analyzed later by a computer with more processing power than the one on 
board a UAV. The second task is essential to the mission and thus is the first priority of 
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the UAV. Each UAV systematically takes data ensuring it receives a data point every 
time it passes around its orbit. The bearing and range values are averaged and this final 
product is uploaded to the gossip network for comparison by the other UAV. Once all 
the UAV have uploaded their data, the orbit phase is complete and the UAV swarm 
moves on according to their mission requirements. 

Once the target’s position is fixed and recorded by the other UAV through the 
gossip network, the UAV swarm then moves to the next target. The second search 
pattern of the UAV swarm is a saw tooth course to the northeast to detect the second 
target and fix its position. This search and seek pattern can be duplicated until the UAV 
swarm reaches its fuel limits or its memory limits with the LPI signal parameters and the 
target’s positions. 



Figure 93. UAV target detection path showing communications with other UAV. 
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VII. CONCLUSIONS 


The four LPI signals introduced in this document perform well when 
implemented in a networked LPI radar system. They illustrate the value of employing 
complex signal designs to improve the performance of LPI radar from three aspects: 
radar capable of identifying and locating targets, the ability to remain undetected during 
periods of radiation and to provide a cooperative signal environment for multiple radar 
sites to function as bistatic radars for one another throughout the network. A networked 
radar system could also use these standard signal design patterns and have another code 
set prepared in reserve for use during times of war. This would reduce the opportunities 
for non-cooperative ESR operators to prosecute the reserve code sets and only present the 
standard design for training and during peacetime. 

Current non-cooperative signal prosecution techniques though capable of 
identifying the presence of a signal are not capable of extracting detailed information of 
regarding the transmitted signal’s parameters. Basic patterns are discernable such as 
approximate center frequency and cycle repetition (for poly-phase only). When 
compared to other poly-phase coded signals (Frank or P1-P4), it is more difficult to 
extract parameter information from the OPS codes. All the signal parameter extraction 
techniques are required to identify these signal parameters. Although these signals do not 
necessarily need to be used by a network of radars, they each have excellent performance 
when used alone in netted LPI radar systems. Longer code sequences could be generated 
to further prevent exploitation by non-cooperative ESR. 

The UAV swarm concept was explored and proved to be a great asset for future 
ISR missions. The capabilities it provides in signal detection and target geo-location 
coupled with removing the risk to personnel has the potential to aid ground forces in 
combat. The UAV swarm can also be used for other missions to include stand-in radar 
jamming and strike missions. The key to the success of the UAV swarm is coordination 
and communication. One UAV is quite capable by itself, but in the same way as the ant 
and termite colonies, the individual is far surpassed by the swarm in a coordinated effort. 
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